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PSEUDOLINEARITY AND EFFICIENCY 
VIA 
DINI DERIVATIVES 
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Department of Mathematics, Miranda House, University of Delhi, Delhi 110007 


AND 


DAVINDER BHATIA 
Department of Mathematics, S.G.T,B. Khalsa College, University of Delhi, Delhi 110007 


(Received 17 May 1988; after revision 15 February 1989; accepted 27 July 1989) 


In this paper, charaterization of pseudolinear functions in terms of Dini 
derivatives is given. Necessary and sufficient conditions for efficiency in terms 
of Dini derivatives are derived for multiobjective programming problems in- 
volving pseudolinear functions. 


INTRODUCTION 


Nonlinear multiobjective programming problems involving pseudolinear functions 
have been studied by Chew and Chool. They have characterized pseudolinear functions 
by means of proportional functions under the assumption of differentiability. 
Rockafellar> has pointed out that the functions involved may not be always differenti- 
able and so characterization of pseudolinear functions without assuming differentia- 
bility is required. 


Kaul et al.4 have defined semilocally pseudolinear functions in terms of right 
derivatives generalizing pseudolinear functions and have extended results of Chew and 
Choo! concerning efficiency for multiobjective programming problems involving 
functions that are semilocally pseudolinear. It was shown be Kaul et al.4 that a function 
fis semilocally pseudolinear on a set [ € A® iff J a real valued function p called pro- 
portional function of f, defined on T x [such that p(x, y) > 0 and f(y) = f(x) 
+ p(x, y) df* (x, y—x) for and x, yin I. 

The purpose of this paper is to define pseudolinear functions in terms of Dini 


derivatives?’ which generalize the class of semilocally pseudolinear functions. In this 
paper, we consider the multiobjective pseudolinear programming problem of the form 


(P) Max f (x) = (f1 (x),..., fe (x)) 


subject to 
eet x= {x€ S | ga (x) = 0ij= PR 
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where it is assumed that ft, = 1,..., k and gj,/ = 1,..., m are real pseudolinear fun- 
ctions defined on a convex subset S of R*. 


Here an alternative characterization of pseudolinearity is provided which makes 
use of Dini Derivatives instead of right derivatives. The useful feature about Dini 
derivatives is that they always exists whereas right derivatives need not. Such a charac- 
terization seems necessary because of the existence of functions of the type given in 


the example presented in section 2, which motivated the authors to make a present 
study. 


Theorems 2 and 3 develop necessary and sufficient conditions for the existence of 
an efficient solution for the above mentioned problem (P), whereas Kaul ef al.4 have 
derived only sufficient condition. 


2. Dint DgrivATIVES AND PSEUDOLINEAR FUNCTIONS 


Definition |—Let f be a real valued function defined over S, a convex subset of 


RN. Letx € RY, v € RV with y7y = 1. Dini derivatives of f in the direction v at 
x are defined as follows : 


De f@) = lim ae {et en fe) | 0 < tm < I/n } 
ae lines 
U : , C+ yr 
D* f (x) = Lowe ont fees esent See mar (= fa I/n } 


Dy“ fix) = lim sup | eae Se) 


-70o {fay —tn 


i eet Pe In} 


Dif (x)= lim inf yeaa Le) 


n>ox {ty} 


O<meiin}. 


u ° 
Here D* f (x) is the upper right derivative, p*' f(x) is the lower right deri- 


. ~U 4 . . as 
vative, D, f (x) is the upper left derivative and D~" f(x) is lower left derivative evalua- 


ted at x in the direction v. Limits can be infinite in the aoe definition. It may easily 
be proved, by using the definitions, that 
(1) Dini derivatives always exist (finite or infinite) for any function f and satisfy 


D*" fix) > D*" f(x), D>” f (x) > D=' f (x) 


Dy (-N@)=— DY s@), DY (—- @) = ~ bt” £@ 
(I If pt f(x) = Dts f(x) (or if D-' f(x) = Dy "f(x)) then the common 


+ 
value, written D » A(x) (or D., f(x)) is just the right (or left) derivative of f at x in the 
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direction y. In general, D. f(x) < Dy f(x) and if D- f(x) = D* f(x), then f 
has the derivative at x in the direction vy and the common value will be denoted by 


Dof (x). 


If f is a function of one variable, then usually, we take the directional vector y 
to be the scalar | and we write 


Di" f (x) = Dt f(x), D¥* f (x) = D*' f(x) 


Dy f(x) = D* f(x), Dy f (x) = Df (x). 
Diewert? introduces pseudoconcave functions as : 


Definition 2—A function f defined over a convex subset S of RY is said to be 
pseudoconcave over S iff for every x° € S,v € R® satisfying 
vy = 1,¢>0,x°+1vE S, Dt" f(x°) < 0° 
implies 
fe +1) <f(). 


In the light of the above definition 2, we may now define pseudoconvex functions 
as follows : 


Definition 3—A function f defined over a convex subset Sof R is said to be 
pseudoconvex over S iff —f is pseudoconcave over S i.e. for every x° € S, v © RN 
satisfying 


Ty = 14> 0, x°+14y ES, Dt f(x°) > 0 implies f (x° + tv) > f(x). 


Hence a pseudolinear function in terms of Dini derivatatives may be defined 
as follows 


Definition 4—A function f defined over a convex subset S of R% is said to be 
pseudolinear over S iff it is both pseudoconcave and pseudoconvex according to defini- 
tions 2 and 3 above. 


Example—Let 


1] , 28x g-— 1 
jor=} 0 eae a, 
1j2atl , 1/2871 < x < 1/2%, n = 0, I, 2, ... 


be a function defined over [—2, If. 
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Some of the Dini derivatives computed at different points of the domain of the 
function are as follows : 


Dt f (— 1) = —e, Dt f(—1) = — 
Du f(0) = 1, Dt! f (0) = 1/2. 
It may be easily verified that 
De f(— 1) <0 = f(—1 +4) <f(- 1),t > 0,(— 1+ 1) €] —1, 9] 
Dt f (0) >0>f(0 +t) > f(0),t > 0,0+ 4 €)0, If. 


Thus fis pseudolinear according to definition 4 over[— 2, I[. It is obvious 
that f does not have the right derivative at x = 0 and hence/ is not semilocally pseu- 
dolinear Kaul et al.4 over [— 2, 1]. 


Assumption—In the sequel, we assume all the functions and Dini derivatives to 
be finite. 


Theorem 1—Let f bea function defined over a convex subset S of RY. Then the 
following statements are equivalent. 


(i) f is pseudolinear over S. 


(ii) There exist real functions p and g defined over S x S such that 
p (x°, x° + tv) > 0, g (x°, x° + tv) > Oand 


f (x? + tv) = f (x°) + p(x°, x° + tv) De ite) 


+ q (x°, x° + 1») Dt £ (x°) (1) 
for any x° € S, v € RN satisfying v’v = 1,t >0,x°+ Wwe S. 


Proor : (i) > (ii). Let x° € S,v € RY satisfying v’v = 1,¢>0,x°+WweEeS 
and f be pseudolinear over S. Therefore 


if D" f (x°) < O then f (x° + tv) < f(x’) 
and if D*’ f (x°) > 0 then f (x° + tv) > f(x) } GZ) 


As the right derivative of f may not exist at every point of S, we have two cases 
Case 1: When right derivative at x° does not exist, i. e. 


u ° u ° 
D*” f (x°) & D*™ F(x?) 
Case 2: When right derivative at x° exists i. e. 


Dy f(x) = Ds" F(x) = D* 4x9, 
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Hence in case (1), relation (2) leads to 
(a) Dt" f(x°) <0 = D* F(x*) < Oand f(x° + tv) < f(x”) 
(b) Dt" f(x) =0 > D* f(x°) < O and f (x° + 1) < f(x’) 
(c) Dt f(x°)> 0 > Dt" f(x°) > Oand f(x° + ty) > f (x°) 
(d) Dt'f(x°) = 0 =D*" f(x°) > 0 and f (x° + t») > f (x°) 


: l 
since D f(x) < 2 f (x) for allx € S and moreover for any x € S, v € RN 
satisfying vv = 1,1>0,x+0ES if f(x +1) =f (x) then D*" f(x) 
ee 
<p if (x)= 0. 


Now we establish statement (ii) in each of the above four possibilities. In pos- 
sibilities (a) and (c) above we may define 


p(x°, x° + tv) = Ble tS SS (x) 
Dt" f (x°) 
and 
q(x, xX += i . Jie Aol 
D, f(x’) 


and the results follows. 


In case of possibility (b), we can define p (x°, x° + tv) to be any positive real 
number and 


Fe) Te) 
Dy’ f (x°) 
Similarly, in case of possibility (d), we define 
f (x° + tv) — f(x’) 
D** f (x°) 


g(x, x° + tv) = 


Pie ek at tv) = 


and g (x°, x° + 1) to be any positive real number. Thus in all the four possibilities, 
statement (ii) holds. 


In case 2, relation (1) can be written as 


f (xe + tr) =f (x) + 2? + 0) D* f (x°) 
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where p’(x°, x° + tv) = p (x°, x° + ty) + 4 (x°, x° +4) and can be defined as 


f pf OP He) = LO ie D* f (x°) #0 


| Das) 
p(x, x°+ tv) = 4 


K, a positive real number, if De ft (x°) = 0. 

The positiveness of the function p’ can be proved exactly in the same manner as 
proved for the functions p and q in Case 1. 

(ii) > (i) : Suppose that there are functions p and g defined over S X S such 
that p (x°, x° + tv) > 0, q (x°, x° + tv) > 0 and relation (1) holds. 

Then Def (x )< O and the fact that ps" ia eA < ps f (x°) imyly that 


ft (x°+ty) < f(x°) showing that fis pseudoconcave over S. Also ptf (x°)> 0 and again 


I , 
the fact that D*" f (x°) > De Sf (x°) imply that f(x° + tv) > f(x°) showing that f 
is pseudoconvex over S. 


Hence / is pseudolinear over S. 
Q.E.D. 


3. NECESSARY AND SUFFICIENT CONDITIONS FOR THE EXISTENCE OF EFFICIENT SOLUTION 


Throughout this section ft, i = 1,..., k, gj, 7 = 1,...,m will be real pseudo- 
linear functions according to definition 4 over a convex subset S of R* with positive 


functions pt, gt,i = 1, ..., k and Py s q, , Jj = |,..., m respectively. 
Consider the following multiobjective pseudolinear programming problem: 


(P) Max f(x) = (fi (x), ..., fe (x)) 


subject to 
FEL Ste © S| Bix) 2 Onset, 2.7). 
For a point x in X, we shall denote by / (x), the set of all 7 such that gj (x} = 0 


Definition 5—A vector v € RN satisfying v7v = 1 is called a feasible direction 
for X at x* € X if there exists a 3 > 0 such that x* + we X forall 0 <1 <3. 


The set of all feasible directions at x* € YX is denoted by D (x*) and 
D (x*) = {v € RY: yTy = 1,98 > Osuch that 
x*+ we Xfor0 <t < 3}. 


Definition 6—A point x* € YX is said to be efficient solution of problem (P) if 
there exists no v € RY with v7y = | such that 


Sa (x* + tv) > felx®) Tet ine 
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and fi (x* + tv) > fi (x*) for at least one f 


where 0 <1? < 6 for some 5 


5 > Osuch thatx* + Wwe YX. 


Note 1: Forx*¢€ X, denote 


Py = Pt (x*, x* + t* v), g*e = ge(x*, x* + 1%), i= 1,...,k 
P'y = P'S (X*.x* + 1% ¥), gS = 4) (xt, x* + t%y), fed, ym 
where we can define? 


ae (x — x*), 
a Te — x*)? (x — xp * =" + VE X for 


== [(x — x*)T (x — x*)]1/2 > 0, 
Lemma-—Let x* be a feasible solution for problem (P), then 


p'; pe gs (x*) + q;" oe gi(x*) > 0 for/ E/ (x*) 
and 


pi D*" gi (x*) +g" D* gy (x*) > —00, 7 & Mx*) implies v € D(x*). 
Proor : Let j/ € / (x*) and ify 


Di: peaks ae De 27 (x*) > O for some 
direction y. 


Then gj (x*) = 0 since j € 1 (x*). Suppose that there exists 8; > 


0 such 
that 0 < t < 83 and gj (x* + 47 v) < 0. Then 
‘* + ‘ l 
pi De" gi (x*) + gi* Di gy (x*) 
* = * —0 
=p’ lim sup 8i (x at v) 0 + gq lim inf BT (x + tv) <= 0 
which contradicts p De gi (x*) + qj De gj (x*) > 0, 
Thus there exists some 57 > 0 such that 
g7(x* + tv) > Ofor0 <t < 8,7 € I (x*). +2{3) 


Now let j & /(x*) and p** D*" gj (x*) + g'* D*" gi (x*) > —co for some direc- 


tiony. Since j & /(x*), therefore gj (x*) > 0. Suppose that there exists 5; > 0 
such that O < tj < 53 and gs (x* + tv) << 0. Then 


? u 7 l * : 
py D, gi(s*)+4;° Dy si (x*) =p; lim sup 


(equation continued on p. 1180) 
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* = * , * t galt bad 
(> 


+ — co 
‘ é ‘ l 
which contradicts py Vase gi (x*) + q;" D* gs (x*) > —©9. 


Thus there exists some 5; > 0 such that 
gy(x* + ty) > 0,0<t G5, fj & I (x*) .-(4) 
(3) and (4) imply that v € D (x*). Q.E.D. 


Theorem 2—If x* is an efficient solution of problem (P) and constraints are 
assumed to satisfy conditions of Lemma then there exists 11 > 0, i =1,...,k and 
pi > 0, j € I (x*) such that 


k 
Su bt Dee) + at DY AO) 


i=] 


‘ ‘ 1 
os > vs [pi* Ds" gs (x*) +. g(* DY gs (x*)]} = 0. ...(5) 
jE1(x*) 


Proor : Let x* be an efficient solution of problem (P) and constraints satisfy the 
conditions of Lemma. Then by Lemma, v € D(x*) and thus there isaé > 0 such 
thatx* + we Xfor0 <t <8, 


We now assert that for 1 <r & k the system 


| jE 1(x*) 6) 

(A) [pr D*" fi (x*) v + 7. D* fi (x*) v]? (x* + tv — x*) > 0, 
| is) sek tser a4 oH 

Lt DE Gty + of DM GEN Gt + — x) >0 
pata) 


has no solution x* + tv € Sfor0<1t< 58. 


Let, if possible, x* + tv € S for some t € ] 0, 8] be a solution of the system A. 
Now as fi’s are pseudolinear over S with positive functions pts and qi’s, therefore 
relation (1) yields 


fa (x* + tv) — fo (x*) = pe (xt, x* + 1) Dt" fi (x*) 


+ qe (x*, x*+ tv) DY fe(x*)i=1,..,k (9) 
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and clearly x = x* + tvE X, 
As in Note 1, we may now define 
th SOS BOS tie 8 
(C Smo WN Caaf 


x = x* + t* vfor7* = [(x — x*)? (x — x*)j1/22 > O, 


Hence (9) implies 
Gla) =— Pp, De fix) + q% DY fi), t= 1, ..,k 


ptt oe. Kae (using (7)): ==3(10) 
Similarly f; (x) — fr (x*) > 0 (using (8)). dL) 


(10) and (11) contradict that x* is efficient for problem (P). Hence system A has 
no solution x* + tv € S forO0 <t<5. Therefore, by Farkas’ Lemma given by 
Mangasarian®, there exist Ar, > 0, wr, 2 0, r = 1, ..., k such that 


Bri [p\* De gs (x*)v + 9‘* pt" gi (x*) vy]? 
JEN(x*) 
k 
a l 
ie S dey (® DY" fet) v + gf De fees) oP 
i=] 
ixr 
— [pt Dt fr(xt)v + af De fe(x*)v? =0. «..(12) 


Summing (12) over r = 1, 2, ..., kK, we get 


k 
> Me [p; Des fa (x*) v + ape D+ty fa (x*)v}F 


i=1 


DS oslo De sscey taf De Gry =o (13) 


jEM(x*) 
where 
k 


At = 1+ Peat 0,. p= 2 pry > O. 
rFi r=1 


Relation (13) yields relation (5). Q.E.D. 


Theorem 3—Let x* be a feasible solution of the problem (P) and let there exist 
Ap 006 1 ey pie 0,7. © I (x*) such that relation (5) holds. Then x* is 


efficient for problem (P). 
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Proor: Let, if possible, x* be not an efficient solution of problem (P). Then 
there exists some v © A% with yTy = 1 such that 
fe (x® + tv) 2 fe (x®), i= 1, -- ,k,iAr 
fr (x* + tv) > fr (X*) 
83 ie 3 ty) Zz 0, J a re ooro m ...(14) 


where 0 < ¢ < 3 for some 3 > 0 such that x* + tv € S. 


s 1% 1% . ee 
Further on using relation (1) and defining Pd, oe q,” as in Note 1, re 


lation (14) gives 


k 
SD ace? Ds FG) + gf DI HEM > o 
ing 
and 
pi Dt" gs (x*) + 4;° Di gs (x*) > 0, J E 1") 
Hence 


: l 
0< SS [pi* D** gy (x*) + a* Dy 84 (x*)] 
jEl(x*) 
k 
‘ l . 
= = Sapt oe per tae ve Ae 
i=1 
<Q 
which is an obvious contradiction. 
Hence x* is efficient for problem (P). 
Q.E.D. 


More general results where in problem (P) the objective functions are 


pseudoconcave and constraints are quasiconcave using Dini derivatives are developed 
in authors subsequent paper. 


Note 2: 
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Let T be a non-empty subset of Zt X Zt where Z+* is the set of positive 
integers and y : T-+>Z+ be a mapping such that for eachn € Zt ,v(xy) =n 
has a finite number of solutions. If f and g are two arithmetic functions, 
then the binary operation } on the set of arithmetic function Fis defined by 


(fi_g(n= > f(x) gly) 
Y (z, Y)=n 


foranyn & Zt. If (F, +, ¥) isa commutative ring, where + denotes the 
usual pointwise addition and % (x, y) # max {x, y} for all (x, y) © T, we 
prove that the ring (F, +, ¥) possesses the unity if and only if for each 
k € Z+, v(x, k) = k has a solution. In such a case the unity can be 
explicitly determined. 


INTRODUCTION 


An arithmetic function is a complex-valued function whose domain is the set of 
positive integers Z*. Let F denote the set of all arithmetic functions. Let T be a 
non-empty subset of Z* x Zt and let 4: T— Z+ be a mapping satisfying the 
following postulates: 


(I) Foreachn € Z*, J (x, y) = nhas a finite number of solutions. 
(II) If (x, y) © 7, then (y, x) € T and v (x, y) = v (y, x). 


(iI) (y,z)€ T and (x, %(y,z))€ T if and only if ‘(x, y) € T and 
(lL (x, y), 2) © T”; whenever one of these conditions holds, we have 


b (x, v (¥, z)) = v (b (2, ¥), 2). 
(IV) (i, 1) = land for each k € Zt, (x, k) = k has a solution. 


(V) Foreach k € Z+, k = max {x € Z+: (x, y) =k for some ye Z*} 
or equivalently (x, y) > max {x, y} for all (x, y) € T. 


If we define the binary operation | on F by 
(f¥yain)= = f(x) g(y) 
y(@-v)=n 
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for eachn € Z*, f, g € F, then using the postulates (I), (II) and (III) it is easily 
seen that (F, +, /) is a commutative ring. 


The following is proved in Sita Ramaiah! : 


Lemma }.1 (Sita Ramaiah!, Lemma 2.1)—An arithmetic function g is an 
identity with respect to if and only if, for any fixed k,n € Zt 


hk f l, ifn=k 
~ &(x)= < PA EG 
y(t-k)=n | O ifaxk. ae 
Let 
Se = {x: (x, k) = k} a tlie) 
ik = min Sz ales) 
and 
== {ik : k = 1, 2, ...}. ...(1.4) 


The number ix in (1.3) exists by postulate (IV). 


If the equation ¥ (x, k) = n has a unique solution in S if n = k and no solution 
in S ifn =k, it was mentioned in Theorem 2 of (Sita Ramaiah!) that the characteristic 
function X of S would be the unity of the commutative ring (F, +, ). 


It was proved in Sita Ramaiah? (Theorem 3.1!) that if g is the unity of (F, +, ¥), 
then g is integer-valued and g (ix) = 1, where ixis as given in (1.3). In addition, if 
g is non-negative, then g must be the characteristic function of the set S_ given in (1.4) 
(See Sita Ramaiah?, Theorem 3.2), Also, an example of aring (/, +, ¥) was given 
(Sita Ramaiah?, Example 3.1) in which the unity assumed negative-values also. 


Suppose g is the unity of (F, +, 4). We shall now investigate the question of 
determining g from the relation (1.1). From (1.1), we have 


a g(x)=1 
xES;, 


so that fork © Sk 


g(k) = 1 — 5 oe (x) ea bee). 
x<k 
xE Sz 

where Sx is as given in (1.2). Let k & Sx. Let xy be the largest element in Sx so that 
Xr < k. From/(l.1), we obtain 


3. gtx) =0. Sate) 
y (ax, )—k 


Now, xr © Sk implies that | (k, xr) = k. Hence from (1.6), 
g(k)=—-— & g (x). sth) 


U(x,Xr) =k 
x<k 
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If g (1) = 1 and g (x) has been defined for 1 < x < k, then the relations (1.5) 
and (1.7) completely determine the value of g(k) forany k € Z*. The task is to 
prove the otherway. We show that (See Lemma 2 2) the relation (1.7) implies that 
g(k) = 0 whenever k & Sx. Using this and (1.5) we ultimately show that (see 
Theorem 2.1) g is the unity of (F, +, /), thus establishing the existence of the unity 
in the ring (F, +, ¥), in the presence of the postulates (I) through (V). 


2. MAIN ResuLts 
First we prove the following : 
Lemma 2.1—We have 
(i) a,b € Sk > v(a, b) E Sx. 
(ii) If a € Sk, then Sa C Sx. 
(iii) (a,b) = k implies that Sa C Sr and Sy C Sx. 


(iv) St = Sxy, where xy is the largest element in Sx. 
Proor : (i) kK = ¢ (a, k) = v (a, U (b, k)) = ¥ (w (a, bd), k). 
(ii) If x € Sa, then $ (x, a) = a. Hence 
¥ (x, k) = 4 (x, 0, k) = (0 (x, a), k) = Va, k) = k, 
so that x € Sx. | | 
(ili) Let b(a,b) =k. Let x © Sa. Then & (x,a) =a. Hence 
¥ (x, k) = ¥ (x, b (a, b)) = ¥( (x, a), b) = Ua, b) = k. 
Hence x € Sz so that Sa C Sx. Similarly Sp C Sx. 


(iv) If k © Sx, then x, = k. So we may assume that k © Sr. Hence 
Xr < k. Since xr € Sz, by (ii), Sxy C Sze. Let x © Sx. Then 
Y (x, xr) € Sk by (i). Also, ) (x, Xr) > xr. Since xy is the largest 


element in Sx, (x, xr) = x, Hence x € Sxy, implying that 
Skt © Sxy. : 


Lemma 2.2—If g is the unity of (F, +, 4), then g (k) = O whenever k & Sk. 


, Proor : If kK € Sx for every k, then there is nothing to prove. We assume that 
& Sh for some k. Let ¢ be the least positive integer such that t & S; clearly 
f>2sincel€ S}. Forl< J< t,7 © Ss. We have from (1.7), 
g(t)= - = g (x) 


Ux,xr=t 
x<t 


where xr is the largest in Sy. 


Now x < ¢ impli ee 
lemma 2.1, ¢ (x, xr) impues that x € Sz. Also, by (iii) of 


=? implies that Ss C S;. Hence x € S:. Since xr € Sz, 
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(i) of Lemma 2.1 implies that t = (x, x7) € St and this is false. Hence the sum on 
the right-hand side defining g(t) is an empty sum. Hence g(t) = 0. We assume 
that g (x) = 0 whenever x & Sz and t < x att 


Let k & Sz. By (1.7), we have 
s(ky)=— g (x) 


(xx J=k 
X= 
where xr is the largest element in Sx. If 1 Sx <tand t (x, xr) = k, then 
x € Sz € Sk, implying that k = & (x, xr) © Sk. We may assume that x > tf. For 
t<x<k, if x & Sz, g (x) = 0, by our induction hypothesis, Arguing as before 
we obtain 
g(k)= — > g (x) = empty sum = 0 
WUx.xr=k 


tSx<k 
xESz 


since k & Sz. We now prove 


Theorem 2.1—Let g be defined by 


f l— f= g(x), if k € Sk 
| x<k 
g (k) = } xES;, 


(2.1) 
| 0, if k & St. 


Then g is the unity of (F, +, ¥). 
Proor : We shall prove that g satisfies (1.1). 


Let xr be the largest element in Sx. By (iv) of Lemma 2.1, St = Sxr. Since 
Xr © Sxr, from (2.1), it is clear that 
l= FS g(x) = 3S g(x). Be ee} 


xESXr xES;, 
It remains to prove that if m and k are positive integers with n < k, then 


> gy) = 0. (23) 
y(v’n)=k 
We distinguish the following cases (in what follows we tacitly assume the results 
of Lemma 2.1): 


Case 1—Let n € Sn. Since g (y) = 0 if y & Sy, we may assume that y € Sy 
in the sum on the left-hand side of (2.3). Now, )(y,) =k implies that Sy C Se 
and Sn C Sx. Since y € Sy and n € Sn, y,n € Sz. Hence k = v(y,n) € Sx. 
Thus we have n € Sn C Skt and k € se. Let Sk = {x1, x2, .. , Xr} with x1 < Xe 


<..<x,=k. Since n © Sn € Sk, we assume that mn = xt where xt € Sz,. 
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First we show that D>» g(y) = 0. Since g(y) = 0, if » & Sy, in this 
Uy, xi) =Xi41 f ; h 
sum, we may assume that y € Sy. Now (y, xt) = Xt41 implies that 


y € Sy © Sz,4, and Sz, ears Since xt € Sz,, we have xt € Sz,,, 5 this 
together with y € Sz,,, implies that x41 = b (y, xt) © Sz,,5- 


Let y © Sz,,, andy & Sz, Since x € Sx,,,, b(y, xt) € Sz,,,. Also, 


u(y, xi) > xt. Hence vb (y, xt) = xt or Xt41, Since xt41 Is the largest element in 
Sz,,, and no element of Sz,,, can exist in the interval (xi, Xi41) as S244 C Sk. 


Hence & (y, xt) = Xi41 since Y(), x1) = x implies that y © Sz,. Therefore by 


1 


(2.2), since Sz, € Sz,,,, we have 


De a ae ee & (y). 
yESz,,, UP X= Xi 41 
V& Sz, 
Let us assume that 
D> g(y)=0 


UY Xi) =Xi+s 
for all swith 1} < s <t<Q&r —i-—J. Weshall prove that 
ose g(y)=0, if xe € Sz, : 
Since 4 (y, x4) = xi, implies that Sz, C Sz,,, 
and xi € Sz, , we have xt © Sasa ‘ 


Ify& Sz, ye Sz,4, then ¢ (y, xt) € Szi4, and %(y, xt) > x4. 


Hence ¥ (y, xt) = xix, or. Therefore 


Oiert 2 AEAT) = S(y) +... + p> g (y). 
yESz;,, yESz;. , yESz;,, 
VESo, WUy,xX )=xXi41 Uy, x)= x14 
...(2.4) 


If xt41 & Sx4,> then b (y, xt) = x¢41 does not occur. Hence the firstsum on the 


right-hand side of (2.4) vanishes. Similar remark applies to the other sums on the 


right-hand side of (2.4). Hence we may assume that xt4s € Sz,,,, for s = 1,2, Dn & 


So, Sri45 C Sz,,,, fors = 1,2, o.0% The variable y in each of the sums on the 
right-hand side of (2.4) can be assumed to be in Sy since g(y) = O if y & Sy. This 
implies that ‘y € Sz,,, need not be mentioned in each of the sums on the Herein 
side of (2.4). Hence (2.4) can be written as 
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o= =f giy= x g(vy) + .0+ x 
yESz;, , Uy Xi )=Xt41 ; YY X= Xi4t1 2 0) 
ye Sz, 
he x g& (y). 
L(y.x¢) = xt 44 

Pad pe -a) 
By our induction hypothesis, the first f — 1 sums on the right side of (2.5) vanish and 
then we obtain 0 = D>) g(y). The induction is complete. Therefore, 


Vy, x= Xe 44 


we have x g (y) = 0, since k = xy. 
Uy. =k 


Case 2—n & Sn. Suppose the equation J(y,m) = k has a solution y with 
y € Sy. Then y € Sy € Sk and Sn € Sk. Let t be the largest element in Sn so 
thatt<n. If ¥(y,n) = k, then we also have 


k — vQ, n) pe vb (y, v(t, n)) zs vy (y, t), n) 


so that & (y, f) is also a solution of the equation » (x,n) =k. Let Yj, Yo, ..., Yr be 
all the elements of Sx which satisfy » (yt. 1) = k, and yi >t for i = 1,2, ..., rand 
yy < Yo <... < Yr. (We may note here that yr = xr the largest element in Sz). 
Also, no ys = ¢t. Forif ye =¢t, then k = ¢ (i,n) = U(t,n) =n. Butn<k. So, 
yt >t for i= 1,2,...,r. We have 


9D IN TF SE TR ay 21 0.) Mie aa a > ). 
y(uon)=k ¥(y.n)=k e(yn)=* 0 y(v'n)=k g() 
¥(vt)=", y(V°t)=Vo y(u't)my, 


Now, U(y, t) = yt and / (jt, 2) = k imply that 
=i) =d(00,9,9) =90,4(40) = 40,9) 
since tf © Sn. Hence we have 


x g@= g(iy+..+ & 2(y) ...(2.6) 
y(y’n) =k ¥ (Vb) ¥y p(ust)=y, 


since t < yt, i = 1, 2, ..., 7 and ¢ © St, each sum on the right-hand side of (2.6) is 
a sum considered in Case | and hence vanishes. Thus 


= gy) =0. 
4 (von)=k 


The proof of Theorem 2.1 is complete. 


We state without proof the following : 


Theorem 2.2—Suppose ¥ satisfies the postulates (1), (11) and (III) of §1 so that 
(F, +,¥) is a commutative ring. If 4(x, y) > max {x, y} for all (x, y) € T, then 
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the unity of (F, +, ) exists if and only if for each k € Zt, equation (x, k) =k 
has a solution. In such a case, the unity is given by (2.1). 


Remark 2.1 : If v (x, y) < max {x, y} for some (x, y) € TJ, then the conclusion 
of Theorem 2.2 need not hold. For example, let T = {(1, 2), (2, 1)} U {(k, k): k 2 2} 
and | on T be defined by | (1, 2) = ¥ (2, 1) = 1 and | (k, k) = k for k > 2. Then ¢ 
satisfies the postulates J, II and III of §1 and clearly for each k € Z*, U(x, k) =k 
has a solution. Note that ) (2,1) = 1 < 2 = max {2,1}. It can be easily shown 
that (for example using (1.1)) that (F, +, 4) does not possess the unity. 


Remark 2.2 : Let T = {(2k, 2k), (2k — 1, 2k), (2k, 2k — 1);k € Zt}. We 
define ): T > Z* by ¢ (x, y) = min {x, y} for all (x, y)E T. It can be shown 
that / satisfies the postulates I, II and III of §1 so that (F, +, ¥) is a commutative 
ring. Also using (1.1) it is not difficult to show that the function of defined by 
g(2k)=1 and g(2k — 1)=0 for k = 1,2,3,..., is the unity of (F, +, ¥). 
Clearly the condition ¢ (x, y) > max {x, y} for all (x, y) € T is violated. 
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The present work has been conceived out of a need to extend the method of 
Ahn to give an iterative procedure for approximating solution of a “quasi- 
linear complementarity problem’ (QLCP). We also give the bound for spec- 
tral radius of them odified matrix in the context of QLCP. Further we give 
more results on the modification of the algorithm of Pang for finding the 
solution of QLCP fora pair (M,q), where M is symmetric and _ positive 
definite. The fixed parameter approach of Pang also has been modified to 
incorporate the variable parameter method in successive iteration process. 


1. INTRODUCTION 


The numerical method for LCP carries with it two methods, e.g, direct and 
indirect methods. Because of the complexity, the use of direct method is restricted 
for large size problems. Therefore iterative methods are well suited for such problems. 
The present work attempts to develop the procedure for finding approximate solutions 
of Quasi-Linear-Complementarity problems by iterative technique. Essentially this ex- 
tends the earlier algorithm of Ahn! used for solution method for Linear Complemen- 
tarity problems. Secondly, we also give an extension of a method of Pang® to 
incorporate QLCP for the same purpose. Our attempt next would be to briefly indicate 
the essential procedure of Pang® because we would refer that latter in our extensions. 


Consider the symmetric LCP (9, M): 
g+Mx>0, x > 0, and x’ (q+ Mx) =0 


where gq © R® and M € R®*" are givenand x € R®. Let(B, C)bea Q-splitting of 
the matrix M, i.e. M = B + Cis a Q-matrix [the LCP (g, B) has a solution for all 


vectors g]. Let E* be a non-negative diagonal matrix with Ey < 1. Define the point 
to set algorithmic map A* as follows: for all vectors x, 

Ak (x*) = solution set of the LCP (q + Cx, B, E* x*). andy ¥ 
The latter LCP (r, B, s) is to find y so that 

r+ By>0, y>sand(y —s)' (r + By) = 0. 
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The LCP (r, B,s) can be converted into the LCP (r + Bs, B) if we translate the 
variable x = y — s since Bis a Q-matrix, the set A* (x*) is non-empty for all vectors 
x. Moreover a vector x* solves the LCP (g, M) if and only if it is fixed point of the 


map A*i.e. x* € A¥ (x*). 

We de fine en iterative technique for solving the LCP (gq, M) given the diagonal 
matrix E* and the Q-splitting (B, C) of the matrix M. Let x® > 0 be an arbitrary non- 
negative vector. In general a* > 0, k > O let x**1 be any vector in the set A*(x*). 

The motivation for using the map A*(x*) lies in the fact that the matrix E* may 
satisfy the bound EN < | after the iteration proceeds onwards after a fixed index ko. 
This idea is compatible with the usual approach in the contraction mapping case where 
a fixed power of a mapping may bea contraction although the original map may not 
be a contraction. 

If Bis a P-matrix (A real matrix A € R”"" is said to bea P-matrix if it has 
positive principal minors) then the set A*(x*) is singleton for all x. In this case, each 
x*+1 will be uniquely defined. 


Pang?! has given necessary and sufficient conditions on the matrix M on the con- 
vergence property, i.e. for all vector q and all starting vector x° > 0, each 
sequence {x*} generated by the iterative technique will converge to some solution of 
the LCP (q, M). 


A Quasi-linear Complementarity (QLCP) can be stated as follows 
Find z € R®” such that 
z— Q0z>0, Mz +q>0(z — Qz)" (Mz + q) = 0. FACT 2) 
With splitting as indicated previously the point to-set algorithmic map takes the shape. 
A (x*) = Solution set of LCP (¢ + BOx*t +C x, 
B, E¥ (1 — Q)x*) (1.3) 
Variantly x* solves the QLCP (1.2) if and only if (1 —.Q) x* isa point in range of 


the set valued map A, i.e. a point in A (x*). The result connected with the conver- 


gence of various dual iterative techniques for the solution of strictly convex quadratic 
program 


min x) =qt T 
ene f(x) =q? x + 4x7 Mx ia) 


can be derived by the methods of Pang®. 
We explain some matrix notations as follows: If A isann x m matrix, « and B 


are subsets of {1,..., n} and il, ..., Mm}, respectively, by Axs we denote the sub-matrix 
of A whose rows and columns are indexed by a and B respectively, Ifa = {1, ..., n}, 
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we denote by Ag the submatrix whose columns of 4 are indexed by 8; similar defini- 
tion applies to Ax. 


2. PRELIMINARIES 
We restate again the QLCP as: 


find z € R”, such that 
z—Qz 20, Mz+qp.0, (z — Qz)’ (Mz + q) = 0 sa(2) 


where Mis ann x nreal and non-symmetric matrix, gism x 1 vector. If we take 
Q = 0 in QLCP we get the same LCP as in Ahn!. 


First of all we describe the notations which occur inthe QLCP. All matrices 
and vectorsare real. A matrix A with m-rows, n-columns is denoted by R”™*”. Row 
i of matrix A is denoted by At and column / by 4; and the element in row i and column 
j by Atj. The transpose of a matrix is denoted by super script 7, such as the trans- 
pose of the matrix A is given by A’, | A | denotes the matrix obtained from the real 
matrix A € K™x® by replacing each element Azz by its absolute value. 


If x € R”, x4 denotes the vector with elements 
(x,); = max {0, xj}; jf = 1, 2,..., m. 
For any x and y in R”, it can be easily shown that 
(i) (x + Ye SX, + 4, 
jij) x y>x%yeSyr. 


A real matrix A € R®”” is said to be a z-matrix (a P-matrix) if it has non-positive off 
diagonal entries (positive primal minors). 


A square matrix with non-positive cff diagonal elements and wiih a non-negative 
inverse in called an M-matrix. It can be easily shown that a matrix which is both a 
Z-matrix and a P-matrix is an M-matrix (or Minkowski matrix). 


Given any real matrix A € R®*”, we define its comparison matrix 


Ac = (Cis) 
by 
Cu = | Att | 
and 
Cy =m j Ag |, tt it = 1,2,085n 


This definition is due to Verga® 
3. ITERATIVE ALGORITHM 


For solving QLCP (2.1) we describe the general fundamental algorithm. 
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Lemma 3.1—Let M € R®*" and E be any positive diagonal matrix, then, 
z— Qz>0, Mz+q>0, (z —Q2z)" (Mz +q)=0 
=z ={(l—Q)z—w E(Mz + q)},, for all or some w > O. 


Its proof is same as in Mangasarian?. This result can be transformed to a fixed 
point problem for solving the equation z = / (z) 


where 
f (2) = (1 — Q)z — wE (Mz + q)}q. 


This result readily leads to the following general algorithm suggested by Mangasarian?. 
We modify this algorithm at certain steps. 


Algorithm 3,1—Let z° > 0, compute 
z6*1 = Q[(1 — Q) zt — wEk (Mzk + q + K* (1 — Q) (zkt1 — z*))}, 
+(1— A) (1 — Q) zt ...(3.1) 
where 
Re 0 1p Oar pe 


and {E*} and {K*} are bounded sequences of matrices in A”x”, with each E* being a 
positive diagonal matrix satisfying E* > « J, for somea > 0 


where / is the identity matrix. 


For the symmetric case Mangasarian has established convergence criteria of this 
general algorithm. We simplify this algorithm by setting 


A= 1, E*¥ = E, K* = K, for each k. 


Remark : As has been indicated in the conclusion, we can relax the above cri- 
teria for fixing the matrix powers E* and K* as constant matricesto derive certain 
variable parameter algorithm as well. 


Algorithm 3.2—Let 29 > 0, compute, 
zh*1 m= [(1 — Q) 2k — wE (Mz* + g + K(1 — Q) (2th -24))), 
ime be Oa ac 


where » > 0, E is a positive diagonal matrix and K is either strictly upper triangular 
or lower triangular matrix. Convergence properties for non-symmetric situations 
can not be established relying on the descent function of the form 


4x" Mx + q?x 


so the recurssive relation between two successive iterations will be utilized here. 
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4. CONVERGENCE PROPERTIES 


First of all we develop the fundamental recursive inequality for Algorithm 3,2 
which will be the basis of convergence. This inequality is derived from the inequality 
properties of x, and y,. 


Lemma 4.1—The kth and (k + 1) th solutions z* and z**1 satisfy the partial 
ordering recurssive inequality : 


| 241 — 2k | < (1 — wE | Kill — | )-1 | (1 —Q) (+ ER) 
— wEM||z* — zk-1|., 4) 
From this Lemma we can produce a condition for the sequence {z*} of Algorithm 
3.2 to be bounded and have an accumulation point which solves the QLCP (2.1). 
If we put O = 0 in (4.1) we have 
tee Pe (Lo BYK 1 )-1, | — o E(M — K)ijz* — zk-1| 
which is the standard form given by Ahn!. 
Theorem 4.1—Suppose that the given iteration parameter w, FE, K and the under- 
lying matrix M satisfy 
H(Z — wE | Kill — Q@|)-1| (1 — Q) U + wERK) — wEM|)< 1 
w0(4.2) 


where » (.) denotes the spectral radius; then the sequence {z*} of Algorithm 3.2 con- 
verges to a solution z* of QLCP. 


The proof is similar to Ahn. 


Here also we can find the same spectral radius which is established by Ahn!, 
simply by taking Q = 0, in (4.2), viz., 


# (U7 — oE| K|)1|7— oE(M— K)|) <1. 


We shall start the next section in which we modify the same result on the con- 
vergence of iterative methods for the symmetric QLCP. Pang had developed necessary 
and sufficient condition (for a fixed parameter) for the convergence of iterative method 
and for solving each individual LCP. We shall extend the method of Pang to incor- 
porate it for the treatment of QLCP. 


5. NON-DEGENERATE CASB 


We classify our analysis into two cases which depends on the nature of the matrix, 
i. e. either the matrix is non-degenerate or positive semi-definite. Since we know that 
the matrix A/ is non-degenerate if all its principal minors are non-zero and the same 
case for the non-degeneracy in linear complementarity theory? the matrix M is non- 
degenerate if and only if the LCP (q, M) has a finite number of solutions for all 


vectors q. 
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The following theorem is the main result of this section. 


Theorem 5.1—Let M be symmetric and non-degenerate matrix. Let (B, C) be 
regular Q-splitting of the matrix M. Let E* bea non-negative diagonal matrix, with 


Ex < 1,foralli. Then the following statements are equivalent: 


(A) for some vector q and any initial vector x9 > 0, any sequence {x*} satisfying 
(1 — Q) x*+1 © A* (x*)is bounded and thus has at least one accumulation point, 


moreover, any such point solves the QLCP (q, M). 

(B) for some vector g, the quadratic function f(x) = g’ x + 4x" Mx is bounded 
below for (1 — O)x > 0. 

(C) for some vector q and any initial vector x® > 0, any sequence {x*} satisfying 
(1— Q) x*¥-1 © A*(x*) converges to solution of the QLCP (q, M). 


Proof can be given in a line of the arguments given in Pang?. 


6. CONCLUDING REMARKS 


As one of the concluding remarks we would like to point out that in the case of 
quasi-linear complementarity problems the algorithm which was developed in section 
3, for the iterative solution technique can as well be generalized for variable parame- 
ters such as the case when the assumptions E = E* and K = K* are relaxed and we 
take uniformly bounded (by matrix norm) matrices E* and K* in the iterative process 
of the algorithm itself. The variable parameter algorithms are still possible to find 
the fixed point for the set-valued maps A (x*). 
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The aim of this paper is to establish some new discrete inequalities in two 
independent variables which can be used as handy tools in the qualitative 
analysis of a new class of finite difference equations involving two independent 
variables. 


1. INTRODUCTION 


The fundamental role played by the discrete inequalities in the development of 
the theory of finite difference equations and numerical analysis is wellknown. A large 
number Of papers dealing with discrete inequalities and their applications have 
appeared during the last few years, see1-19 and some of the references given therein. 
Although stimulating research works have been undertaken in this direction, there 
are still a number of interesting classes of multidimensional finite difference equations 
which needs new types of discrete inequalities in their analysis. Our objective here is 
to present some new discrete inequalities in two independent variables which can be 
used as handy tools in the qualitative analysis of a new class of finite difference 
equations in two independent variables. In order to convey the importance of our 
results to the literature, we present applications of some of our inequalities to the 
study of boundedness, uniqueness and continuous dependence of the solutions of a 
new class of fourth order finite difference equations in two independent variables. 


2. STATEMENT OF RESULTS 


We first summarise some basic notations and definitions which will be used 
n—1 

throughout this paper. Let No = {0, 1,2, ..}. The expression u (0) + a b (s) 
3= 


represents a solution of the linear difference equation Au (mn) = b(n) for n€ No, 


where A is the operator defined by Au (n) =u(n + 1) — u(n), The expression 
m4 

u(0) !i b(s) represents a solution of the linear difference equation u(n + 1) 
s=0 


= b(n) u(n) forn € No. We use the usual convention of writing oi b(s) = 0 
Ss 
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and II 6b(s) = 1, if ® is the empty set. We also use the the following notations of 
the ReariSts 

Ai u(m,n) = u(m + 1, n) — u (m,n), 

Agu (m,n) = u(m,n + 1) — u(m, n) 
for m,n © No. We often use the letters mandn to denote the two independent 
variables which are members of No. 


For convenience we list the following hypotheses : 


(Hi) u(m,n) and A(m,n) are real-valued nonnegative functions defined for 
m,n € No. 


(He) pi(m, n), po (m,n), p3 (m,n) are real-valued positive functions defined for 
mn € No. 


(Hs) a (m, n) is real-valued, positive and nondecreasing function in both the variables 
m and n in No. 


(H4) u (m,n) > up > 0, up is a constant, h (m,n) > 0 are real-valued functions 
defined for m,n € No. 


(Hs) g(u) is continuous, nondecreasing real-valued function defined on an Nepietoe 
I = [uo, 2°), uo > O is a constant, and g(u) > 0 on (uo, °°), g (uo) = 


(He) 91 (m,n), g2 (m,n), q3 (m,n) are real-valued positive functions defined for 
m,n € No. 


(H7) W (wu) is continuous, nondecreasing and submultiplicative real-valued function 
defined on an interval J, and W (u) > 0 on (uo, e-), W (uo) = 0. 


A useful two independent variable discrete inequality is embodied in the 
following theorem. 


Theorem 1—Suppose (Hi) and (He) are true. If 


m_} z~1 n-l 


eC I 
A tN EES ee Pi (x, n) Sexe ey 7) > wes rs y) 
z= y=0 
v=]. 
x SS h (s, 1) u(s, 1) Al) 
t-0 


for m,n © No, where c is a nonnegative constant, then 


—1 


n-1 
vimose TT [14 ‘ 7 P3 (Ss, y) 


t= 
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y=1 


x h (s, » | wate) 


t=0 
for myn © No. 


A slightly different version of Theorem | is given in the following theorem. 


Theorem 2—Suppose (H1), (Ha) and (Hs) are true. If 


m1 z-1 n-l 
oe A cas a Pi ay n) 1) P3 ea y) 
z=0 6=0 v=0 
y-1 
x h(s, t) u (s, t) 3) 


for m,n € No, then 


n- 


£1 1 
u (m, n) < a (m, n) hee nore “pi (x, 1) Sten) sce y) 
om y~ 


r=0 


y_1 


x h (s, t) ...(4) 
Soy 
for m,n € No. 


Another interesting and useful discrete inequality is established in the following 
theorem. 


Theorem 3—Suppose (Ha), (H4) and (Hs) are true. If 


m-}\ pat n-1 
] ] ] 
BS dr iio > Di (x, n) Sao Sy ps (Ss, ¥) 
x=0 e=0 i= 
y-1 
x h (s, t) g (u (s, t)) (>) 


for m,n € No, where c is a nonnegative constant, then forO < m < m,0 GC n<m, 
m, m4, n,n € No, 
nd 


m=1. @-1 
1 l eee. 
u(m,n) < G-1 | G(@) + Sore p2 (s, n) > ps (s, Y) 
x=0 s=0 yO 


y-1 
se SS h(s, 0 | ...(6) 


t= 
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where 
caf. 

G(r) = [ ay ,r > ug with ro > uo nog 8 

7 


0 


G~! is the inverse of G and mj, 11 € No are chosen so that 


m-1 2-1 n-1 y-1 
] ] 
LC hres Pa 
(oe Py (x, n) P2 (Ss, n) Ps (Ss, y) 
2=0 s=0 y70 tO 


x h(s, t) € Dom (G-1) 
for m,n € No and0<ma<am,0<n< ny. 


We next establish the following more general inequality which may be convenient 
in some applications. 


Theorem 4—Suppose (H:), (He), (He) and (Hz) are true. If 


m1 z-1 n-) 
] l ] 
u(m,n)<et+ S at a 5 tall ma > ‘a wee 
“5 Pi (x, n) P2 (Ss, n) P3 (s, y) 
r=0 s=0 y=0 


y-1 
- SS h (s, t) u(s, t) 
t=0 
m-_-j Gal n-1 
+ Sate’ abe Sato 
qi (x, 7) 92 (s, n) q3 (Ss, y) 
z=0 s™0 v=0 
y-1 
x ay k (s, t) W(u(s, t)) ...(8) 


t=) 


for m,n © No, wherecisa nonnegative constant and k (m, n)is a real-valued non- 
negative function defined for m,n € No, thenforO< me m2, 0 | n & no, m, mo, 
n,n2z © No 


m1 z-1 
u(m,n) < Q(m,n)Q-1! a( ot. pea Eerie Ta 
: L () 2S qu (x, n) > q2 (Ss, n) 
nS ' y-1 
xX ener k (s, nwo (s, t)) | ...(9) 
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where 
m_} Lae | nm-1 
ome = [1+ nord nba D 
| Pr (x, n) P2 (8, n) a3 y) 
xz=0 s=0 v0 
Pa 
x h(s,t)) ...(10) 
t=0 
and 


Tr 
d a ; 
an=-| 7%. r > uo with ro > uo 


ol & | 1) 
fi) 
Q-1 is the inverse of Q and mo, no © No are chosen so that 


z-1 n—1 


10+ Satay n) BED D we BEN 


s=0 


y-1 
k (t, s) W(Q (t, s)) © Dom (9"1) 


t=0 
for msn © No and0O gt m<m,0K€n<n 


3. Proots oF THEOREMS 1—4 
In order to establish the inequality (2) in Theorem |, we first assume that c > 0 


and define a function z (m,n) by 
n—1 


plenies © Ss Rie Pi (x, 1) Saas “ath n) Dat ps (Ss, ¥) 
far 
y-1 
x h (s, t) u (s, t). gs(12) 
From (12) it is easy to observe that 
z (0,2) = z (m, 0) = foto} 
and 
n-} ori 
pi (m, n) Ay z (m,n) = h (s, t) u (s, t) 
ee po is y) <—/ P3 ries 


s=0 y 
... (14) 
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"1 y-1 
1 

po (m, n) Ar [pi (m, n) Aa z (m,n) = >. md; h(m, t) u(m, t) 

é : BILL? 
n-1 

p3 (m, n) Ae [pe (m, n) Ax [pi (m, n) Ay z (m, n)]] = mS h(m, t) u(m, t) 
| ...(16) 
As [ps (m, n) Ae [pe (m, n) Ay [pi (m, n) Ay z (m, n)J]] = A (m, n) u are 


Using the fact that (uv, m,n) & z(m, n) in (17) we have 


As (p3 (m, n) Ae [pe (m, n) Ay [pr (m, n) Ay z (m, n)]]] < A (m, n) z (m, i 
ny, 


From the definition of u (m, n) we observe that z (m,n) <z(m,n + 1) form, n E Np. 
Using this fact in (18) we see that 


Ps (m,n + 1) Ag[p2 (m,n + 1) Aq [p31 (myn +1) At 2 (m, n+)]] 


z(m,n + 1) 
P3(m, n) Ae [pe (m, n) Ay [pi (m, n) Ay z (m, ny] < h(m, n). 
a (m, I 1) as Kis 


From (19) and the fact that ps3 (m,n) Ae [po (m, n) My [py (m, n) Ay z (m, n)]] > 0 
from (16), we observe that 





_P3a(m, n + 1) Ae [pe (m,n + 1) dy [pi (m,n + 1) A, z(m,n + 1)]] 
Z(m,n + 1) 


p3 (m, n) As [po (m, n) Ar [pi (m, n) Ay z (m, n)]] 
z (m, n) 





<A (m, n). 
...(20) 


Now keeping m fixed in (20), set n = t and sum over ¢ = 0,1,2,...,2—1 and use 


the fact that ps (m, 0) Ao [pe (m, 0) Ai [pi (m, 0) Ay z (m, 0)]] = 0, from (16), to 
Obtain the estimate ; 


P3 (m, n) Ae [pe (m, n) Ay [pi (m, n) Ay z(m,n)]] _ 
Pees) als i 21 GeO 
t=0 


...(21) 
From (2!) andin view of the facts that z (m,n) <z(m,n + 1) and P2 (m, n) 
Ai [pi (m, n) Ay z (m, n)] > 0, we observe that 
P2 (m 


a+ 1) Ay [px (m, n + 1) Ay z (m,n + 1)} 
aa ee nt) a 8) 
Z(m,n+ 1) 
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_ P2 (m, n) 44 [pi (m, n) O12 (m, n)) 
z (m, n) 
n-] 
h (m, t). A he} 


MO n) 
1-0 


Keeping m fixed in (22) set n = y and sum over y = 0, 1, 2, ..., n--1 and use the fact 
that po (m, 0) a1 [pi (m, 0) Ay z (m, 0)] = 0 from (15), to obtain the estimate 


mes y- 
P2 (m, n) Ay [pi (m, n) Ay z (m, n) 
z (m, n) “Saar Ps IER y) 5 ged 


(23) 


From (23) and in view of the facts that z(m,n)<z(m-+1,n) and pj; (m,n) 
Ai; z(m, n) > O from (14), we observe that 





“pi(m + 1,1) A, z(m + 1,7) _ pi m,n) Ay z (m, n) 
z(m + E n) z(m, n) 
n-1 y-1 
A gore oe), aa S h (m, 1). 24) 


Now keeping n fixed in (24), set m = s and sum over s = 0, 1,2, .., m — 1 and use 
the fact that p; (0, n) Ai z (0, n) = 0 from (14), to obtain the estimate 


m-] n-1 y-1 
Ay z (m, n) l > aE l >. 
, < ee h (s, t). 
Z (m, n) P1 (m, n) P2 (s, “pe (s, n) P3 aaa r) 
s=0 y=0 t=0 


oF ys), 
From (25) we see that 


m_ 
Rect eds 
z(m + 1,n) &2z(m, ny aay Pi Cae p2(s, n) > Ps (s, ¥) 
ho +53 


y-1 


x yy h(s, n |. -..(26) 


t=0 


Now keeping n fixed in (26), set n= x and substitute x = 0,1,2,..,m-— 1 
successively and use the fact that z (0, n) = ¢ from (13), to obtain the estimate 


n-1 


2-1 
z(m,n)&c It - —o n) ELS ae n) a mE 
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y-1 


x > 2,9 |. 


t=0 


Substituting this bound on z (m, n) on the right side of (1) we obtain the inequality 
in (2). 


Now suppose c = 0. Then from (1) we see that the inequality 


n—1 
u(m,n) < «+ Sata Pl = n) Sain a n) aga math y) 
y-1 
x h (s, t) u (s, t) 


holds for every arbitrary positive number « and m,n € No, which by the above 
argument yields the estimate 


m1 


u(m,n) < ibe te PI a “Py (x, n) Sain ED Sain ea 
SS his.) |. ixlae) 


Since u (m,n) > 0 and € > 0 is arbitrary number independent of m,n then from 
(27) it follows that u (m,n) = 0. This completes the proof of Theorem 1. 


Since a (m, n) is positive and nondecreasing, we observe from (3) that 


nl 

u (m, n) >; ! ] 

a (m, n) mee Pr ace n) 4 — Pz (s, n) Y Ds (Ss, Y) 
v=0 


Sy h (s, t) ri. an ‘ 


t=0 





Now an application of Theorem | yields the required bound in (4) and the proof of 
Theorem 2 is complete. 

In order to establish the inequality (6) in Theorem 3, let « > 0 and u, (m, n) 
= u(m,n) + «> up forall m,n € No. Then from (5) we see that 
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m1 2-1 ni 
ue (m,n) Se +e + Serre cen Se 
Pi (x, n) p2(s, n) P3 (s, y) 

2=0 s=0 yv=0 


y-1 
x SS h (s, t) g (us (s, t) — €) 
t=0 
m_j 
> Te at yy Pi = n) aw Sake p2(s, n) yor =n y) 
y-1 
x h (s, t) & (ue (Ss, t)). .. (28) 
t=0 
Define a function z (m, n) by 
m-1 @-1 ' n-1 , 
sey aati Saray 2 pr (x, m2) Po (s, n) Day (5,9) 
Te 
x h (s, t) g (ue (8, 1). (29) 
t=0 
From (29) it is easy to observe that 
z(m,0)=z(0,n)=ct+e ...(30) 
and 
m4 n-1 , y-1 
pi (m, n) Ay z (m,n) = Dp 2 Pe > h (8, #)8 (use(s, £)). 
pto.) 
nl 
pa (m,n) &:[p1 (m,n) d1 2 (m,n) = > new Km 1) 8 (ue (m,t)) 
v=0 
oH f4) 


n-1 
ps (m,n) Az [p2 (m,n) O1[p1 (m,n) O1 2 (m, n)]] = real a, Bea (m1) 
FAEKS: 


Ao[p3 (m,n) Ae [po (m,n) Ar [ps (m,n) M1 z (m, n)]l] = h(n, m) g (ue (m, 0). 
...(34) 
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Using the fact that u. (™, n) < z (m,n) in (34) we have 


Ao [pz (m,n) Se [po (m,n) Ay [pi (m, n) Ay z (m, n)]]] < h(m, n) g (z (m, n)). 
35) 


From the definition of z(m,7n) in (29) we observe that z (m,n) < z (m,n + 1) for 
m,n € No. Using this and the fact that 


ps3 (m, n) Ae [po (m, n) Ay [pi (m, n) Ay z (m, n)}] > 0 
from (33), we observe from (35) that 


_P3 (m,n + 1) Ag [po (m,n + 1) Ay [pi (m,n + 1) Adz (m,n + IDI 
&(z(m,n + 1)) 





__ Pa (m,n) Ae [pe (m, n) Sy [pi (m, n) Ay z (m, n))] 
g (z (m, n)) 





< h(m, n). 
...(36) 


Now by following exactly the same steps as in the proof of Theorem | below the 
inequality (20) up to the inequality (25) with suitable changes, we obtain 


nj y-1 


4iz (m,n) — 
g (Zz (m, T) aoe Pl on n) xa > Po i n) > Fa P3  (s, ‘y) rm h(s, t). 
y—- 








Pes ed 
From (7) and (37) we have 
2(m41,N) d A ( ) 
G (z (m + 1,n)) — G(z(m, n)) = y 1 Z (m,n) 
(m + ) (z (m, n)) me WES < Bee paps 
m-1 n=t 
PN ESS es, 
px (m, n) p2(s, n) Ps (s, y) 
&=0 v=0 
y-1 
x h (s, t). Be 
t=0 


Now keeping n fixed in (38), set m — x and sum over x = 0, 1, 2,,....m — 1 to 
obtain the estimate 


m-1 ny n-1 


G > G € ; 
(z (m,n) < G(e + 6) + Dron) n) " p2 (s, n) Sian 
pm 


s=0 


de 


x dS h (s, ). .-.(39) 


t=0 
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The bound in (6) now follows by substituting the bound for z (m,n) from (39) in (28) 


and letting « + 0. The subintervals of No for m, n are obvious and the proof of 
Theorem 3 is complete. 


In order to prove the inequality (9) in Theorem 4, lete > 0 and wu, (m,n) 
=u(m,n) + € > up form,n € No. Then from (8) we see that 


t-1 n-1 


ue (m,n) < tet Say Tea n) rey y) 


v=0 


< 
=" 


x 
7M; 
o 


m_1 z-1 


h(s, t) (ue (s, t) ae) int 





=r m) “qa oo n) 


z=0 
n_1 y-1 
l 
x we k (s, t) W (u (s, t) — 6) 
y=0 t=0 
m_1 2-1 n- 
eS SS 
aa hs Pi (x, n) P2 (Ss, n) ps (s, y) 
z=0 s=0 yv=0 
y-1 
x Ah (3.2) ucts; t) + Ser ae mi Seam say “ 
t=0 x=0 
ni y-1 
c Ss. EetV WE (es. 2): ...(40) 
P3 (s, y) 
y=0 t=0 
Define 
z-1 n-1 
BAmtgB) im Cots €1r Sato Fae “1 (x, Nn) . POEL Re we AOE y) 
a 
y-1 
x x k(s, t) W (us (s, t)) (41) 
t=0 
then (40) can be restated as 
m1 G=1 I 
l 3 l 
ue (m,n) & a (m,n) + > Pi a Pi (x, n) Ls P2(s,n) ps (s, Y) 
2=0 8=0 y-0 


a 
x5 h(s, t) Us (s, t). 
t=0 
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Since a (m, n) is positive and nondecreasing function in both the variables m and n, we 
have from Theorem 2 

ue (m,n) < a(m,n) O(m, n) ...(42) 
where Q (m, n) is as defined in (10). Since W is submultiplicative, we have 

W (us (m, n)) < W (a(m, n)) W (O (m, n)). ...(43) 
From (41) and (43) we have 


w-1 


m- 1 al 1 
1 l >) ie ap 
aimn<ctet S —O oo Sai q3 (S, ¥) 
z=0 8-0 =? 


y-1 ) 
x & k(s, t)W (QO (s, t)) W (a (s, t)). 
t=0 


Now by following the proof of Theorem 3 with suitable modifications we obtain 
m-1 r-1 ? 
] 
- Q-l Setitgenn sho) obs cae 
2 (on) OG ise) Sere g2 (s, 0) 
xz=0 8=0 
n—] 


y-1 
] 
S ae s W , .. (44 


yv=0 


The desired bound in (9)now follows by substituting (44) in (42) and letting « > 0. 


The subintervals of No for m and n are obvious. This completes the proof of Theorem 
4. 


4. Some APPLICATIONS 


In this section, we present some applications of our results to the study of boun- 
dedness, uniqueness and continuous dependence of the solutions of a new class of non- 
linear finite difference equations in two independent variables. Fach of these applica- 


tions could be stated formally asa theorem. This has not been done so as not to ob- 
scure the essential ideas with technical details. 


Example 1—As a first application, we obtain a bound on the solution of a 
nonlinear fourth order finite difference equation 


Az [az (m, n) As [ao (m, n) Ay [a1 (m,n) Aq u (m, n)\]] = f (m, n, u) (m, n)) 


osx (49) 
with the given boundary conditions at m = 0,n=0 


u (0, n) = ¢4 (n) 
a1 (0,n) Ay u (0, n) = $2 (n) 
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az (m, 0) Ay [a1 (m, 0) Ay u (m, 0)] = di (m) 
az (m, 0) Ag [az (m, 0) Ay [a1 (m, 0) Ay u (m, 0)]] = ve (m). ...(46) 


ape 3 9 9 
Here a1, a2, a3 are real valued positive functions defined on N) ,f:N) X R > R, 


where R denotes the set of real numbers; ¢1 (7), ¢2 (”), Yi (m), ve (m) are real-valued 
nonnegative functions defined for m,n € No. We assume that 


| f (m,n, u) | SAC, n) | u | ...(47) 


where /' (m, n) is a real-valued nonnegative function defined for m, n € No. Itis 
easy to observe that the problem (45) — (46) is equivalent to the equation 


m-1 G1 n— 
1 
pe er) > ay (x, n) Dee a ag (Ss, n) : Ge y) 
“=O s=0 y-0 
y-1 
x'S f(s, tu (s, 1) (48) 
t=0 


where 


sat m-1 
' 
b (m,n) = $1 (n) + SS pny ie Aas Seen 
z=) z=0 


cr | m-1 £-1 
] l l 
ss Patent: is) by ay (x, n) > ag (x, n) 
e=0 z=0 s=0 
n- 
x te (s) =, hay i (49) 
heed § 


Suppose that 
|b(mn)| Sk ...(50) 


where k is a nonnegative constant. Using (47), (50) in (48) we have 
n-1 


m_1 z-1 | 
] ee ee 
|u(m,n)| Sk + Say D> sor az (s, Y) 
xz=0 80 Ho 


rie | 
mS en (set) te Ge, tb): 
t=0 


Now an application of Theorem | yields the bound on the solution u (m, n) of (45)-(46) 


in terms of the known functions. 
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Example 2—Asa second application, we shall discuss the uniqueness of the 
solution of the problem (45)— (46). We assume that the function / in (45) satisfies 


| f (m,n, u) — f (m,n, i) | Shum n)|u—a| Pe 


where / (m, n) is asin Example |. The problem (45)—(46) is equivalent to the equa- 
tion (48). Then for any two solutions uv and a of (45)—(46) we have 


z=] 


m_1 i 
a he Fee, 
|u(m,n) — i(m,n)| Se + -s=—ea D> a2 (s, n) 
z=0 s=0 


n. 


_ 


y-1 
| S . 2 
x as (s, y) ra h (s, t) | u (s, t) oe (s, t) | -a(32) 


—) 


Ve 


where « > O is arbitrary constant. The assumption (51) is used to get the inequality 
in (52). Now an application of Theorem 1 yields 
| u (m,n) — a (m,n) | 


n—1 


n-1 z-1 , y-1 
_ l 
ee ilar Sc rer) roe a kl 


=0 


—~ 


Since « > 0 is arbitrary we have u = di. e. there is at most one solution of the pro- 
blem (45) — (46). 


Example 3—Our third application is an example of continuous dependence of 
the solution on the equation and boundary data. Consider the problem (45)—(46) in 
Example | and the problem 


Ag [ag (mm, n) As [az (m, n) Ay [ay (m, n) Ay z (m,n)}|| = F (m, n, z (m, n)) 


with the given boundary conditions at m = 0,n = 0 Se 
z(0,n) = ¢o (n) . 
a1 (0, n) Ay z (0, n) = yo (n) 
a2 (m, 0) Ay [ay (m, 0) Ay z (m, 0)]] = py (m) 
a3 (m, 0) Az [az (m, 0) Ay [ai (m, 0) Ay z (m, 0)]] = de (mn), ...(54) 


Here 41, ao, ag are as in Example 1, F: Ni x R>R, $1 (n), bo (n), dy (m), be (m) 


are real-valued nonnegative functions defined for m,n € No. The equations equiva- 
lent to (45) — (46) and (53) — (54) are (48) and 
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m i Traik: n-] 
z(m,n) = b(m,n) + > > >= 
z=0 " ( n) 3-0 ig (s, ") hm! sf ts, ) 
y i 
x BF (s, tz (8, 0) ..(55) 


where b (m, n) is obtained from the definition of b (m, n) by replacing $1 (1), by (n), U4 


(m), v2 (m) in the right side in (49) by dy (n), bo (n), v4 (m), vs (m) respectively. From 
(48) and (55) we have 


u (m,n) — z(m, n) 
= b(m,n) — b(m, n) 


m_} z-1 y 


1 
=. | , 
; > ay (Xx, ae ae ca ny > ag (Ss, y) 


s=0 v=0 


-1 
XE (£65, 1,U(s, 9) — Fs, t,2 (5, 0} ...(56) 


Suppose that the function / in (45) satisfies the condition (51) and further we assume 
that 


| b(m,n) —b (m,n) | <e AAT. 
m 1 b steomt | “-1 é y—1 
» ne S aa (s,m) >, GW ut SS, Leela) 

a AS, £2348, 0)) | me ...(58) 


where ¢ > 0 is arbitrary constant. By substracting and adding f (s, ft, z (s, t)) in the 
braces on the right side of equation (56) and using (51), (57), (58), we obtain 





m=) Epa | 
] ] 
lum m—zimnis2+ Sao > wear 
x=0 s=0 
n-1 y-1 
h(s, t) |u(s, t) — z(s, t) |. (oo) 
v=0 t-0 


Now an application of Theorem 1 yields 


| u (m,n) — z(m, n) | 
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a oe n-1 
ele tl E sf resp <i aoe 

y-1 

> : (so | ...(60) 


If h (m, n) is bounded on some compact set 0 < m < mo, 0 <n < no, mM, mo, Nn, No 
€ No, then the quantity in braces on the right in (60) is bounded by some constant 
M on the settO <m €m,0<n <n. Therefore | u (m,n) — z(m,n)| {2M 
on the setO0 < m < mo, 0 Sn € no; so the solution u (m, n) of (45) — (46) depends 
continuously on f and the boundary data. If «0, then|wu (m,n) — z(m,n) | 
— 0 on this set. 
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Results on existence of solutions to Periodic Boundary Value Problems for 
an infinite System of nonlinear second order differential equationshave been 
discussed and a uniqueness result is presented. 


1. INTRODUCTION 


There is a large literature on the existence of solutions for periodic boundary 
value problems (PBVP’s for short) of nonlinear scalar differential equations?7~10°12"14, 
Results on the existence of solutions for first and second order PBVP’s have been 
obtained earlier?’®'1914 by combining two basic techniques, namely the method of 
upper and lower solutions and the alternative method. In deed, a number of existence 
theorems for periodic solutions have been obtained by the technique of upper and lower 
solutions, and others have been and are being obtained by the alternative method. 
Kannan and Lakshmikantham®” showed that the use of a result proved by Cesari and 
Kannan? by the alternative method could remarkably improve the arguments by upper 
and lower solutions. New existence theorems could be obtained and previous 
ones could be obtained by a more uniform approach. Further work in this direc- 
tion was done by Rao and Vatsalal®. For an exposition on the alternative 
method, though preceeding the work in Cesari and Kannan*, we refer to Cesari?. 
Also we refer the readers to Bernfeld and Lakshnikantham! for a general treat- 
ment of boundary value problems. Extremal solutions for nonlinear boundary value 
problems have been obtained by many workers?’5°8’9'°11°16'l7 by employing mono- 
tone iterative scheme and this method is found to be constructive. PBVP’s for infinite 
systems of first and second order equations have been studied earlier$11’1415'17 and 
these results extend the results obtained for scalar equations. 


In this paper we consider, as it was done earlier 45-716, an infinite system of 
second order ordinary differential equations with periodic boundary conditions and 
discuss the questions of existence and uniqueness of solutions again by combining the 
method of upper and lower solutions with the result in Cesari and Kannan® from the 
alternative method. The organization of the present paper is as follows : Our notations 
and terminology are fairly consistent and can be understood by referring to earlier 
workers®’11'15°16, However for the sake of completeness we describe them briefly in 
section 2. Section 3 deals with existence results employing the alternative method. In 
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section 4, we develop a monotone iterative technique to obtain coupled extremal 
quasi-solutions for these systems. The uniqueness results are discussed in section 5. 


2. NOTATIONS AND TERMINOLOGY 


Let J be the interval [0, 27] and E= R® = R x Rx ... in which R = (—°9, 
cc). Also Jet Zt denote the set of all positive integers. 


We consider the PBVP 
=U fas ft (1, Wy ue jy t. CS, SeGuie acatoe le 


u (0) = u (27), u' (0) = wu’ (27) 
where f: / X E x R => E iscontinuous. * 


For each i € Z*, we define two sets Pt and Qi such that Z* — {i} Ps U Qs. Here- 
after P(g) denotes the generic element of Pi (Qi) respectively whenever the set is 
nonempty. Moreover the vectoruw € Emay be written as u = (ut, [u]p, [u]o), for 
pe Pi, gE GQ. Then the PBVP (2.1) becomes 


—ui = fit (t, ui, [u)p, (ula, ui ) epee’ 


u (0) = uw (27), u' (0) = uw’ (27). 


Without further mention we assume thati € Z* and all the inequalities between 
vectors hold component wise. 


We state the following assumptions which will be used in our subsequent dis- 
cussion : 


(Ao) «, 8 € C2[/, E), a(t) B(t), te 7 

(41) (i) « (0) = « (2x), a’ (0) > x’ (Qn) and — a" <f (t,o, a") 
for all o such that « (1) € o < 8 (r) and ot = of (t), 1€ (0, 2n] 
(ii) B (0) = 8 (27), 8’ (0) < B’ (27) and pe > fe (1,9, B, ) 
for all o such that « (t) € « < A(t) and ct = Bs (t), 1€ (0, 27] 

(42) Fort€ sath < u(t) <B(t)andvw € R we have 


he(|ve|) if| re] <a 


lft(t,u, ul < | 
he (a) if [ue] > 


for some dt > e¢ = 1/2m max {at (0) — Be (2m) | ms 
— $ ’ «(2 — 0 : : 
(0, cc) is continuous for each i. | 2¢ (27) — Br (0) |}, At: [0, cc) 
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Also there exists N > 0 depending only on «, 8 and A such that 


Y 

sds . 
| hi (s) > ae Bi (t) ae Sey “4 (t). 
t 


where 
lt = min {di, Nt} 


(As) fis completely continuous on / x E xX R. 


(44) Fort € La(th<v<u<B(t)and | z|< di 
fi (t, ut, [u]p, [ula. zt) — fa (t, vt, [u]p, [ul], 24) > — Me (ut — vi) 


for some Mi: > 0. Where di = max { Ni, max | a (t)|, max |B’ (¢) | 
I I ; 


(As) f possess a mixed quasi-monotone property (mq mp) that is ft (t, us, [u]p, [u]a, zt) 
is monotone nondecreasing in [v]p and monotone nonincreasing in [u]g. 


The functions «, @ € C2 [/, E] with « (t) < 6 (t) on/ are said to be coupled 
lower and upper quasi-solutions of (2.1) respectively if 


— a; <M (1, 24, far, [Ble 2), (0) = + (2m), #0) > &! Cn) 


— 8; > fa(t, Br [Ele (ale, 8; ), 8 0) = B (2n), B (0) < B' (2n). 
The functions x, » € C2[J/, E] are said to be coupled quasi-solutions of (2.1) if 


, ~ hs («, xt, [x]p, [y]a, x; ) , x (0) = x (27), x’ (0) = x’ (27) 


—-y=ft (:. yt, Ly]o, [x]a, 9; ) y (0) = » (27), y' (0) = y’ (27). 


In the special case where all Qi’s are empty, quasi-solutions are just solutions and in 
the case Pi’s are empty, the quasi-solutions that result are most useful since they may 
be obtained most easily. We can also define coupled minimal and maximal quasi-solu- 


tions analogously. 


We state the lemma which is a modified version of a known result. 


Lemma 2.1—Let the assumptions (4o) and (Ag) hold, then for any solution 
u € C2 [I, E] of (2.1) with « (t) < u(t) <& 6 (1) on J, we have 


|u'(t)| < Non/. 


The proof of this lemma follows from Lemma 1.1 of Das and Devasahayam. 
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3, EXISTENCE RESULTS 


In this section we discuss the existence of solutions of PBVP’s by the device, al- 
ready used many times in the method of upper and lower solutions, of defining the 
function F as follows : 


Fi (t, u, vi) = ft (t, p(t, u), vt) + re (t, u) 


where 
Pi (t, u) = max {as (t), min (ui, Bt (t))} 
and 
( LBL Ea if ui > Be (t) 
| 1 + ul 
2 
Bae : 0, if at (t) < um < Be (ft) 
| Nb lee ,ifus < a¢ (t). 
1 + ui; 


Consider the following PBVP 


aoe PR (:, u, U; ), u (0) =u (27), u’ (0) = u' (27). Be PH 
Lemma 3.1—Let (Ao) and (A1) hold and let u be a solution of (3.1). Then 
a(t)< u(t) < B(t) onl. 


Proor : First we claim that « (1) < u(t) on J. Suppose not, then we can find 
ato « Jand ane > O such that for some k € Zt 


ak (fo) = uk (to) +e, at (th u(t)+e tE I. Walore) 


If to € (0, 27), we have «, (to) = uw, (fo) and a’. (to) < ux (to). From (3.2) ax (to) 


> uk (to) and hence px (to, u (to)) = ak (to). In view of (Ay) (i) and using the defini- 
tion of F we have y 


Sk (to 9, a (to)) > — at (to) 


7 


— Uy (to) = Fr (to, u (to), uj, (to) 
> fk (to, P (to u (to)), ao (10). 


Since « (t) < p(t, u (t)) 8 (t) and Pk (to, u (f0)) = 


: ak (fo), ere) 
by choosing « = p (to, u (to)). (to), we get a contradiction 
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If to = 0, from (3.2) we obtain 


ak (0) = uz (0) + € = ak (27) and a, (0) € uj, (0) anda, (27) > uw, (27), 
since «x (0) = ak (2x), wx (0) = ux (27). 


And in view of (A1) (i) it follows that «/ (27) = wu), (27) and as before 
Si (27, 5, «) (27)) = —a) (27) = — uf, (2m) = Fr (2n, u (27), ut, (27) 


> fk (27, p (27, u (27)), a, (27)). 
Since « (27) = p (27, u (27)) < B (27) and px (27, u (27)) = «% (27), we again get a 
contradiction by choosing 6 = p (27, u (27)). On similar lines we can prove that u(t) 


< 8 (ft) on / and this completes the proof. 


Lemma 3.2—Let the assumptions (Ag) — (Ag) hold. Then there exist «0, 29 such 
that the following are true : 


(43 ) bobs © C2IL El, wo (1) < Bott), re I 
(4F ) @ 20.) = 2 Cn), 4% (0) > % (2m) and 
aon Fil t, Ga. ) for « such that 
ao (t) S 5 < Bo(t) and ot = a, (t), t © (0, 2x]. 
(ii) 80 (0) = Bo (27), By (0) & By (27) and 
= F(t, eee ) for all ¢ such that 


ug (t) < ¢ < Bo(t) and ct = Bo, (t), t E (0, 27). 


(43 ) The condition (42) holds with F replacing f and with respect to the pair 
(xo, Bo). 


Proor : Let at > 0, bt > 0 for alli € Z* be numbers and define ao,¢ (t) = «a 
(t) — at and Bo,t (t) = Bi (1) + bt. Then it is easy to see following the proof of Lemma 


2.2 Rao and Vatsalal5 that {| Aj and (4 ) hold. However for the sake of complete- 


ness we establish ( Ay ) ; 


Fort € J, «(t) <u(t)<B(thandu; € R, we have 
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[Fi (t, vw r= 1e(im a) 
m(im 1), if | 3f6) alee a 


he (dt), if |v’, | > de. 


\ 
a aN 


Since hi (s) is a positive constant for s > di, there exists an N r > dt such that 





N* N; 

sds sds 
i Gy | in (8) > ap Bo, (t) — ere ao,t (1) 
sere 


where 


* 


; = min {ei, ai} and a4 = 


c 


= max [ | «o,¢ (0) — Boe (27) |, 


| ao,¢ (2) — Bo,¢ (0) | ]. 
This proves that (42) holds with Ft replacing fi, i € Z*. 


For o, « such that 29 (1) < a, o< Bo (t), ot = ao,¢ and of = Bot, define 








| Fe(t,o, ee ) + Boe ite > Boy 
| l+u; 
| t 
| 
| ww — aoe [ - 
: Fi( 1c R! ee I 
= » Korg , Ort 
Gi (t, u) = * Bot — cost 
+ Fy (+, a, ee if xo,¢ < Bo,t 


; | 
’ aot —-ut. 

| Fi (1, S, %,, ) + ——— if uw < ag,. 

* 


| 1+u3 


Since <o,¢ < Bo,t for all ¢ © J, Gs (t, u) is well defined. We now consider the modi- 
fied PBVP 


— uy = Gt (t, u), uw (0) = u (27), uw’ (0) = wu’ (27). yatsia) 


Lemma 3.3—Assume that (Ao), (A1) and (As) hold. 


: Then the problem (3.3) has 
a unique solution u satisfying 


age (ft) cu(r)< S Bot (t), t € I. 
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Proor : It is clear from the definition that Gi (t, u) is completely continuous and 
bounded on / x E. Hence we can find a positive number J that depends on ag and 
Bo such that ||G (t, u)|| C J. Let ¥ = Lo (0, 2] define Lin = — u; and then D (ZL) 


= {9 € X: 9, 9’ are real valued absolutely continuous on [0, 27], ep EX, 9 (=o 
(27) and 9’ (0) = 9 (27)}. Let OS?’ be the nonlinear operator defined by 


Wu = Gt, u). 
Then the BVP (3.3) may be translated into the operator equation 
Lu= Wu 


Notice that Xo, the Kernel of L consists of all constant functions and hence X1 where 
X = Xo @ Xi is the class of all vector functions whose average on [0, 27) is zero. 
Also we can define the operators P and H satisfying the conditions of Theorem 2.1 
of Kannan and Lakshmikantham’. Since G is bounded we can find constants A and B 
that depend only on 9, 8 such that any solution of (2.2) in Theorem 2.1 of Kannan 


and Lakshmikantham’ satisfies |lw1 (t)|| < A and Ilu; (|| < Bfor allt € J. 
Hence by Theorem 2.1 of Kannan and Lakshmikantham’ it is enough to find an Ro,i 
> 0 such that 

< WD (uot + 1,1), uot > > Oor < O ...(2.4) 
for all | vot | = Rot and |u| < A,|u,, | < Bonsi€ Z*. Since Xo consists 
of all constants functions and uo € Xo, (3.4) can be written as 


Qa 
[ Ge (s, Ro + u1(s))ds <0 
0 


and 


Qr 
{ Gi(s, — Ro + m (s)) ds 20 
0 
choose Ro,i > 0 large enough so that Ro,t + w1,¢ > Max Boi (t) and — Rot + u1,8 


< min «4,¢(t), Using the definition of Gi and ( At )we see that 
I 


‘On 


Qn 
| Gi(s, Ro + um (s)ds < | FES, %; Boss (s)) ds < 0 
0 0 
and 
Qn Qy7 


fo (s, — Ro +m (s))ds > | Ft (s, 4, %,, (s))ds > 0 
0 0 
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for any arbitrary but fixed o, o satisfying 


a9 (t) <6, @ < By (t), ot ag, (t) and ot = Bo,t (t) 


Hence by Theorem 2,1 of Kannan and Lakshmikantham’, there exists a solution u to 
the problem (3.3). Proceeding on the similar lines of that of Lemma 3.1, it is easy to 
see that u(t) satisfies the inequality 


ao,t (t) {us (t) < Boyt (t), ¢ © J. From this it follows that for each i € 


Z*, uw is a solution of 


pee ui — 40,8 = F a ! 
~ Mt “Boe — ost E (. % Por _ (« %> Org )] 
+ Fe (1,9, 04, ) was(3.5) 
u (O) = u (2m), u’ (0) = u’ (27) 


where a (t) < 5, ¢ < Bo (1), ot = a0, (t) and ot = Bo, (t). Following the proof of 
Lemma 2.3 of Rao and Vatsala!5 it can be shown that the solution u(t) of (3.5) is 
unique. 


We now prove our main result of this section. 


Theorem 3,1—Let the assumptions (Ao) — (A3) be satisfied. Then the PBVP (2.1) 
has a solution u such that « (¢) < u(t) < B(t) and | u(t)| < Non I, where N de- 
pends only on «, 8 and the Nagume function. 


Proor : Consider the boundary value problem 


ul = (1.1, u' ), u (0) = u (27), u’ (0) = w’ (27) 136) 
where 


He (1, 0, 4; Jaa ret uu )+0-a G(,0,2 € (0, 1] 


In view of Lemma 3.2, one can verify that xo and Bp satisfy ( Ay ) : ( A= ) and 


( A ) with respect to Hi. Using the arguments similar to that of Lemma 3.], it is 
easy to show that, if v),s for \ € (0, 1) is a solution of (3.6) then 40,¢ (t) & uy,¢ (1), 


S Bos (t), t€ J, ie Zt. Since Hi satisfies a Nagumo condition, we have 
| U4 (t)| =C,t © IJ where the constant C is independent of A. In view of Lemma 3.3, 
it follows that for A € [0,1) all possible solutions of (3.6) satisfy ap (1) & uy (t) Sfo(t) 
and | u. (t)| < C forall t € /. Also forA = 0, the problem (3.6) has a unique 
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solution. Thus one can choose a bounded, closed convex set B in the (u, u') space H1 
[0, 27) = Le [0, 27] such that (3.6) has no solution on the boundary of B for 
AE (0, 1) and it has a unique solution in the interior of B for A = 0. Hence by Leray- 
Schauder theory, the problem (3.6) has a solutionu for A = 1. By Lemma 3.1 we 
have « (t) < u(t) P(t), t © Jand this wv is a solution of (2.1) satisfying « (t) <u 
(t) < B(t), ¢ © FZ. Thus (Ag) implies that | u’ (t) | < N,t © I where Nis the Nagume 
constant vector. This completes the proof. 


As a special case of Theorem 3.1 we have the following result. 


Corollary 3.1—Let the assumptions (Ao), (A2), (As) and (As) hold. Further 
assume that there exist coupled lower and upper quasi-solutions « and §. Then the 
conclusion of Theorem 3.1 is true. 


Remark 3.1 : Note that if x, 8 are coupled lower and upper quasi-solutions, then 
(Ai) holds if f has mg mp and also that (A1) implies «, B are coupled lower and upper 
quasi-solutions of (2.1). 


Our results can also be extended to the following infinite system of second order 
equations with homogeneous Neumann boundary conditions (NBVP for short). That 
is 


<u =fi( uu, ), uO =w Qn) =0 3.7) 


where 
fe Cl[lx EX R, E] andi € Z*. 


Theorem 3.2—Let (Ao), (Ag) and (A3) hold. Further let the following condition 
hold: 


(B) (i) 2% (0) > 0, « (27) < O and 
== Ge. (0). aft (1, a, &, ),t € (0,27] for all « such that 
a(t)<o<BP(t)ands =a(t), tEL71€ ih 
(ii) 6’ (0) < 0, B' (27) > Oand 
—-R ah (+, o, B; Ne t € (0, 2z] for all 
such that a(t) < o < @ (t) and «1 = Bi (1), t€ 1,i€ Zt. 


Then the problem (3.7) has a solution u such that 
a(t)<u(t)sS B(t)and|u'(t)| S N on I. 


Where N depends only on ~, B and the Nagume function. | | 
The proof is similar to the proof of Theorem (3.1) with appropriate modifications. 
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4. MONOTONE ITERATIVE MeTHOD 
For any yn,» © C[/, E], «(t) & y(t), uv (t) < B (tf), t © J, we consider the 
quasilinear PBVP 
— ub = Gt (t, ut, [u]p, [u]a, uf), u (0) = u (27), wu’ (0) = u' (27) ...(4.1) 
where 
Ge (t, ut, [uJp, [uJo, uy ) = fe (t, nt, [y]p, [u]a, g (uy ))—Me (us — xs)) 


and 


g (u, ) = max | — a, min (u; a) | 


Notice that G is defined on / x [«, B] x R and (A4) is equivalent to 
( At )f (1, ut, [u]p, (ula, g (u; )) = felt, vt []p, [u]e, 
g ( x; )) > — Mi (uw — ve) 


for Mi>0,r € J, a(t}< vu < 6 (t) and uw E R, 


Relative to the PBVP (4.1) we prove the following lemmas. 


Lemma 4,1—Let the assumptions (Ao) — (Ag), (Aa) and (As) be satisfied. Then 
the assumptions (41) and (Ao) are true with respect to the PVBP (4.1). Thatis «, 8 


are also coupled lower and upper quasi-solutions of the PBVP (4.1) and Gi satisfies the 
modified Nagumo condition relative to @. By 


Proor : Using the arguments of Lemma 3.1 of Lakshmikantham et al.1! it is 


easy to show that «, B are also coupled lower and upper quasi solutions of the PBVP 
(4.1). However when (42) holds we have 


me ( Iw | Jit tu ica 
| fe fF ut, [u]p, [ula, us Lite { 
{ he (a)if tul | > de ... (4.2) 


fort € /andsomed > g@ = = max { | o¢ (0) — Bs (2m) | , | at (2m) — Be (0) ] } 


*(1)Su(t)< 6 (t),u © Rand © C [[0, co), (0, ©0)), also 
\ 


sds 
ares > rahe st 
% 
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where 

kk = min {ds, M4}. 
From this it follows that any solution wu € C2 [/, E] of (2.1) satisfies lias) eh Se Ne 
te J. 


For? € /, a(t)<u(t) < B(t) (onus (t) and u; € R we have 


| Mi(u — 1) |< Mev 


where 

y= eR Bi (t) — Fai ag (ft) 
Further 

| Ge (t, ut [u], [u]o, ui ) | = He Clu"). 
where 


. h M. 5 ‘=o 
Ki (s) = a tyifs <a 
hi(di) + Meyifs > de 


Evidently since Hi (s) is a positive constant for s > dt, there exists an ee > di such 
that 








N* N° 
i i 
{ s ds = | Sis.” = : 
Hi(s) ~ Ai (s) ~ t 
&% q% 


and this proves that Gi also satisfies (Az). Hence the proof is complete. 

We now prove a result on existence and uniqueness of solutions of the PBVP 
(4.1). 

Lemma 4.2—Let the assumptions (Ao) — (As) hold. Then there exists a solution 
u of (4.1) such that a(t) G u(t) (B (t)and | uw; (t) | < Mi on /. Further more the 
solution u (tf) is unique. 


Proor : By Lemma (4.1) we have all the assumptions (Ag) — (Ag) satisfied with 
respect to the PBVP (4.1). Hence by Theorem (3.1), there exists a solution u(t) of 


(4.1) with « (t) G u(t) < B(t) and |u; (t)| < Neon J, Using the arguments similar 
to that of Lemma 3.2 of Lakshmikantham et al.11 we can show that the solution wu (ft) is 


unique. 
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Since for every y, » € [«, B], the PBVP (4.1) has a unique solution u, we define 
the mapping A by 


A(y, ») =u .--(4.3) 
and study the properties of this mapping in the next lemma. 


Lemma 4.3—Under the assumptions of Lemma 4,2, the mapping A defined by 
(4.3) satisfies the following properties. 


(i) « < A (a, B)and 8 > A (8, a) 
(ii) ForaC yc p< B. A(a, 2) SA (u, 9). 


PRoor : We shall only prove that B > A (f, «) since similar arguments prove that 
a & A (a, §), 
Let A (B, 2) = wu, where uw is the unique solution of the PBVP (4.1) with » = B and 
# = 4% Let p(s) = u(t) — B(t). Suppose that the inequality p(t) <0, t€ / is 
false. Then there exists a 19 € / and ane > O such that for some index k € Zt, we 
have 


pk (fo) = « and g (t) <e forallt€ Jandie Zt. ...(4.4) 


If fo € (0, 27), we have 9) (to) = 0 and ?;. (to) < 0. 
Also 


& (uy, (to)) = g(8z (to)) = Bf (to). 
At ¢ = 10 using (A,) and (4.1), we have 


0 > 9% (to) = ut (to) — 8* (to) 


> — Gk (to, uk (to), [4 (to)]p, [4 (to)]a,u?. (t0)) 


+ Se (to, Bx (to) [B (to)}p, [x (to), Bi, (to)) 
= Mke > 0, a contradiction. 
If to = 0, then using the boundary conditions, (4) (ii) and (4.4) we obtain 
Pk (0) = ux (0) — Be (0) = ux (2m) — Be (27) = pk (27) = 
?, (0) € 0 and 9, (27) > 0. 
Also 


9, (0) = wi, (0) — B, (0) > uy, (2m) — B, (27) = 9 (2z). 
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Hence 9, (27) = 0 and using (A) and (4.1) we get 


Q, (27) = ui. (27) — By (27) > Mk € > 0 which is again a contradiction. 


To prove (ii) let 4, 7 € [«, 8) such thatyn Sp. Let A (1, uw) = x, A(u,y) = y and 
v(t) =x(t)— y(t). Ifthe inequality | (t) <O fort € Jis false, then there exist 
fo € Jand ane > 0 such that for some k € Z*, we have 


ok (to) = « and tt (t) Ce fort € Jandi € Z*. ...(4.5) 
If to € (0, 277), we have 


wi (to) = 0 and vy (to) < 0. 
From (4.5), (4; and (As) and inview of the definition of G, 
0 2 ¥ (to) = x, (to) — yz (to) 
= — fk (to, nk, (np, Lula, 8 (x), (to)) + Me (xe (to) —nk (to)) 


+ fk (to, pk, [u]p, (nla, & (¥;, (to)) — Mk (ve (to— px (to) 
= Mk (xk (to) — ye (to) > 0, a contradiction. 
If to = 0, then from (4.5) and the boundary conditions we obtain vx (27) =e and 
vi, (27) = 0 and as before we get a contradiction at ¢ = 27. This completes the proof 
The following is the main theorem of this section. 


Theorem 4.1—Let the assumptions (Ag) — (A5) be satisfied. Then there exist 
monotone sequences {un (t)}, {8n(t)} with 79 = «, Bo = B such that on (t) and Bn (ft) 
converge uniformly and monotonically to Pp (t) and r (1) respectively on J. Also (P, r) 
are coupled minimal and maximal quasi-solutions of the PBVP (2.1). More precisely, 
if (x, y) are any coupled quasi-solutions of (2.1) satisfying « < x, y < 8, then 


i=ojS a aee &... Sons... QPaxyyars... S Ba 
< ... © Bi < Bo = 8 on I. -«-(4,6) 


Further more any other solution u of (2.1) satisfying « (t) < u < B (1) on / also satis- 
fies (4.6) on J. 

Proor : We know from Lemma 4.2 that for any %, » € [2, §], the PBVP (4.1) 
has a unique solution u (t) such thata(t)<uSB (t) and|]uy (t) | S Mon J, 


A . 
where N; is the Nagumo constant relative to G. In view of Lemma 4.3 we may define 
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the sequences an = A (%n-1, Bn-1) and Bn = A (Bn 1, 7n-1), 1 = |, 2, 3, ... such that 


ay = a and Bp) = 8 and an < Bn foreachn. Since «9 < 71<f1 < fo, by induction 
and the arguments similar to those used in Lemma 4.3, we can establish that {an}, {Bn) 


are monotone sequences such that 
19 0, <... Can < Bn < ... B2 < Bi < Bo on /. 
Where an (t) and Bn (1) satisfy 


— an, = ft (t, an-rt, [an-r]p, [Bn-1)a, g (a, )) 


— Mt (ant — yn-1,1), an (0) = on (27), 0° (0) = % (2m) 
pl4ed) 


— Bog = ft (ts Bait, [Bn-a)p, [on-rake, g ( Pn )) 


—Mt (Bn,t — Bn1-1), Bn (0) = Bn (27), BY, (0) = Bi (2m) ...(4.8) 
and 
LMI.) BOL<™. 
From (4.7), 


t 0 
— ont (t) = J J J4 (8, amast (5), lama]p, [2n-1]a, 8 (ah, (5))) ds do 


- 
oa 


! Mt (an, (8) — any, (s)) ds do + Cnt + Ang. 


Since for each i, fi 1s completely continuous, the sequence {«n} is uniformly bounded 
and equi-continuous. Thus {a,} contains a subsequence which is uniformly conv ; 

by the Arzela-Ascoli theorem. In view of the fact that {an} is Ae.» deers 
Sequence ccnverges uniformly on /. Further more the uniform boundedness - es: 


sequence {«” } implies that the sequence {a’ } is equi-continuous and uniformly bou- 


nded. Thus {«) } contains a subsequence which is uniformly convergent on /. Similar 


reasoning for {Bn}. Thus we can find sub-sequences which we again denote by {an} 
nS, 


{8n} converging uniformly and monotonically to P,r. Where P (t) = lim « (7) 
. - 
andr(t)= lim Bn(t), t€ J. And this implies that ee 
no 


— Pi fa (t, Pa, [P]p, [r]a, g (P5, )), (0) = P (27), P’ (0) = p’ (27) 


~ ri = feb re Ur, [Pe g (rf )), (0) = 7 (2m), ¢ (0) =r’ (2) 
(4.9) 
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Following a continuation argument similar to that of [Bernfeld and Lakshmtkantham! 

p. 32], one can prove that f and r are actually coupled quasi-solutions of th PBVP (2:1), 

If (u,v) are any coupled solutions of (2.1) such u,v € [x, B]t © Zand | u’ |, |v’ | 
A 


<N <d on/, employing induction principle and the arguments similar to those used 
earlier, it can be shown that #n < u,v < 8non/. Hence we have P <u,v Cron /, 
proving that P, r are coupled minimal and maximal quasi-solutions of the PBVP (2.1). 


Since any solution u of (2.1) satisfying « < u< B on / may be regarded as (u, u) 
coupled quasi-solution of (2.1), we also have P <u Sr on /. This completes the proof 
of the theorem. 


5. UNIQUENESS RFSULT 
We now present a result on the uniqueness of solutions of the PBVP (2.1), 
Theorem 5.1—Assume (Ao) — (A3). In addition for each i € Z*, there exists a 
constant Z; > 0 such that 


(uc — ve) [ft (t, uu, uy )— fe (tv, ¥, < — Le (us — v4)? Berk bei 8 


whenever « (t) < u,v < 2(t)t © Jand uy, — = 0. Then the PBVP (2.1) has a 
unique solution u (ft) satisfying « (t)< u(t) < B(t) on /. 


Proor : By Theorem 3.! we know that the PBVP (2.1) has a solution. If pos- 
sible, let u and v be two Solutions for the PBVP (2.1) satisfying «(t)< u (t), v(t) 
< B6(t)fort € /. 


We define 

pt (t) = (us (t) — ve (t))? 
and observe that p+ (0) = pt (27), H (0) = p; (27) for eachi € Z* 
and 


ni (t) = —2(ue() —u (OLA yw I-fi(ay % 


+2 (4 (1) —¥, (OP tS) 


We claim that ps (t) = 0 for all i © Z* on / If not there exists a to € J and ane > 0 
such that for some k € Z* 


uk (to) = cand pt (t) < forallt € 7, i € Z*. cad bee 


If to € (0, 27), we have 


py (to) = Oand ut S 0. 
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Thus from (5.2) and (5.1) we have 


O > pe (to) = — 2 (uk (to) — ve (to)) [ fe (to, u, uw, (to)) 


iv Sk (to, Vy, Vi (to))] nig 2 (u, (to) =, “s (to))? 
=> 2 Lk pk (to) > 0, a contradiction. 


If to = 0, we obtain pe (0) = ¢ = pk (27) and He (0) < Oand pw, (27) > 0.How- 


ever p), (0) = p», (27) and hence», (0) = pw, (27) =O and consequently »; (A) 
< 0 for A = 0, 2r. 
Also 


O > py (A) = 2Lk pe (A) > 0, for A = 0, 27 


which is again a contradiction. Hence the proof is complete. 


Corollary 5.1—Assuming the conditions of Theorem 4.1 and the hypothesis (5.1) 
one may conclude the existence of a unique solution for the PBVP (2.1). In this case, 
if suffices to show that the coupled minimal and maximal quasi-solutions p (t) and r (ft) 
are identical. This can be accomplished by defining the function 


ut (t) = (Pa (t) — re (t))? for eachi € Zt andr € J 


and proceeding along the lines of the proof of Theorem 5.1. 
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The connection between the 3—/ and the 6—/ coefficients to a set of six F (1) 
sand a set of three (or, equivalently a set of four) nS (1) s, respectively, 


is used to obtain sets of Regge 3 x 3 and Bargmann—Shelepin 4 x 3 
symbels, Closed form expressicns are obtained for the polynomial zeros 
of degree n of these coefficients. 


INTRODUCTION 


In literature, classical symmetries of the 3—j and the 6—j coefficients are 12 
and 24 in number. By relating the 3—/ coefficient toa 3 x 3 magic square symbol, 
Regge! showed that it has 72 symmetries. Bargmann? and Shelepin® related the 6—/ 
coefficient to a 4 x 3 symbol, which exhibits all the 144 symmetries (including the 
classical symmetries) discovered by Regge’. We have shown that sets of p41Fp (1) s 
are necessary and sufficient to account for all the known symmetries of the 3—j and 
the 6—j coefficients—explicitly, while a set of six 3Fo (1) s represent the 3-j co- 
efficient®, either a set | of three or an equivalent set II of four 4F3 (1) s represent the 
6 —-j coefficient®’. Here we establish a connection between these sets of generalized 
hypergeometric functions and sets of Regge, Bargmann-Shelepin symbols for the 3-j 
and the 6—/ coefficients, respectively. Using these, closed form expressions have been 
obtained for the polyncmial zeros of degree n of the 3—j and 6—j coefficients, which 


have been the subject of detailed Study-especially when. = 1 or 2—by several authors 
in recent years. 


CLOSED Form ExprEssiONs 

The 3—j coefficient has been defined by Wigner® as : 
ji je J 44-jo7m 3 

( = 8(m + me + ms3,0)(—1)' *? 3A jejs) {(it 
my, m2 m3 a | 


+ ma)! (jt — mi)N2 x DT (— 1)? (21 0 {z — ak)! 
3 : “1 
I (B:—z)!} 1 ost t} 
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where 
max (a1, %2) < z < min (61, 82, Bs) 
“1 = fi — Ja + me, a2 = jo — fg — my 
Bi = j1 — m, Bo = je + mo, B3 = Ji + jo — js. ii(2) 
and 


A(xyz)=((—x+y4+z)!(e —yt 2)l(x + py — 2)!/ 
(x +ytz+ ipl hie) 


and 5 (x, y) is the Kronecker delta function. Regge! discovered new symmetries by 
associating the 3 — / coefficient with a magic 3 x 3 square symbol : 


fhAti2 tis A-—jetis titse - js 


I te San 
es me mg ) seth feed | my J2 — me St oie 
ji +m J2 + me Ja + ms 
= [Reel (4) 


and asserted that the 3 x 3 square symbol represents the invariance of the 3 —/ coeffi- 
cient to 3! column and 3! row permutations and toa reflection about its diagonal. 
Thus, Regge! established the existence of a 72-element symmetry group, comprising 
the well-known classical symmetries (which arise due to column permutations and to 
the space reflection : mi > — mtarising dueto the interchange of the second and 
third rows of (4)) and six new symmetries known as Regge symmetries of the 3—/ 
coefficient. 


It has been shown by one of us (Srinivasa Rao®) that the 3—/ coefficient can 
be represented by a set of six 3F2 (1) s: 


hi jz Js 3 
( ) = 3 (my + me + m3, 0)(— leer) oT (Ralf + Dy? 
my, m2 m3 c bokeh 


x (f (1 — A, 1 — B, 1 — C, D, E)}"} 3F2 (A, B, C; D, E; 1) 


(5) 
where 
A = — Rop, B= — Rza, C = — Rir, D = Rar —Rop + 1, 
= Ror ar R3a = 1; he (x, y; — = I (x) ry) core 
J=fitse + Js ---(6) 
and 
Rsp — Roa for even permutations 
3 (pgr) = { 
R3p — Rea + J for odd permutations 
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for all permutations of (p g r) = (123). Since each one of the six 3F2 (1) 5 represents 
only 12 symmetries arising from its invariance to 3! numerator and 2! denominator 
parameter permutations, the set of six 3F2 (1) s is necessary and sufficient to account 
for all the 72 symmetries of the 3—/ coefficient. 


Using the properties of the elements of the 3 x 3 square symbol : 
Rip a Rmp = Rna ao Rnr ret ie 4 


for (/mn) and ( pqr) being (123) cyclically, and the defining relations (6) for the nume- 
rator and denominator parameters, one can easily show that: 


SB Dit tl gee dees iC ae 
| 

Reel) = ety =: CdD She oe 
Ri Fiky 8 Ag aD Ate 


From (8) it is straightforward to obtain the closed form expression : 


ieee js ai agi ree (— B—C+ D —1)/2 
(A-C+E—1)/22 (~B+C-—D+1)/2 


. My Me Ms 


(- A—B+ D+ E—2)/2 9 
(- A+ B+ D — E)2 S 
The parameters A and B are negative integers by definition. So, if welet P= — A 
and 0 = — Band set C = — 1, we get, for the polynomial zeros of degree one the 
expression : 
( hh je js ) Ny ( (P + E)/2 (Q + D)/2 
my m2 m3 (=P E)/2 8 DIZ 
(P+Q+D+E- 2)/2 
-Ata0 
(P—Q@+ D— E)/2 ) (1) 
with the constraint equation : 
PQ = DE : ofl) 


which is a homogeneous multiplicative Diophantine equation (see, Bell9) of degree 2 


whose complete solutions can be expressed in terms of four parameters. Identifying 
the four parameter solution to be : 


P =ab,Q = cd, D=bdandE=ac yay be) 


and substituting it in (10), we obtain the result given by Brudno!l9, The polynomial 
zeros of degree n arise due to the truncation of the 3F¢ (1) series (5). By setting anyone 
of the numerator parameters, (say C) to — n and equating the sum of the (n 5 1) 
terms to zero, one obtains the constraint equation which must be satisfied by the 
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numerator and denominator parameters of the 3Fo (1) for realizing the polynomial 
zeros of degree n. 


The 6—/j coefficient has been expressed by Regge‘ to be: 


iy 


N 2 (— 1)? (p + 1)! tl “pt 
= ee ! ae y ; Ee : 
def a CPSU NY TP esa) HL (Bs — py? (13) 


with 
N=A(abe)A(cde)A(acf) A(bdf) 
ma =-atb+ea=-c+dt+eo=atectfou=b+dt+f 
Bi=a+b+c+d, Ppo=at+dt+et+f, Pp=btctet+f 


and max (a1, %2, «3, %4) < p & min (Bj, Bo, 63). The 144 symmetries of the 6—/ 
coefficient arise due to the invariance of (13) to the 4! permutations of the «’s and the 
3! permutations of the {’s. This is explicit in the notation of Bargmann? and Shelepin® 
wherein : 
Pri ay 2 Pe = xy 1 6a — 103 
alee By — &e Be — ae Bs — ag 


hep ie = || Rex] (14) 


Bi — 73 Bo—ag B3 —- ag 





Bi — ag Bo — 44 «3 — a 


which is invariant to 4! row permutations and 3! column permutations. The elements 
of ||Rix|| satisfy the 18 relations : 


Rez + Rmn = Rkn + Rmk 
and 

Rak + Rmn = Ran +Rmk (15) 
fork = mand k snandk, morn being 1, 2 or 3. Equivalently, every 2 x 2 co- 


factor in (14), sav | i, : | satisfies the condition: « + 5 = B+ y. 


We® have shown that the 6-—-/ coefficient can be represented by a set / of three 
4aF3(1)s as: 


me == (—1)2t1NT (| — E) (TU —4,1—8, 1—C, 1 — D, 


dcf 
Rays als (A, B,; C; D; E, FG, 1) mA) 
where 
A= — Rip, B = — Rop, C= — Rap, D= — Rap, 
E = — Rip—Rop — Roc — Rar — 1, F = Roe — Rap + 1, G 


= Rar — Rap + 1 +(17) 
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for (p qr) = (123) cyclic; and use has been made of (15) in arriving at (17). It is now 
possible to express the standard Bargmann-Shelepin 4 x 3 symbol in terms of the 
numerator and denominator parameters of the set / of 4F3 (1) s, using (15) again, as: 


ey is ee es i 
[abe apt p/p ayia = oa 
} Ae Loh Fue Co 1G Dee 
4)D)-FisD dlaGreteent 


= || Rix ...(18) 


where the negative denominator parameter E does not appear in (18) and the 4 x 3 
symbol in (18) exhibits only 48 symmetries which arise due to the invariance of the 
6—j coefficient to 4! row permutations and 2! column permutations in ||Rex||. The 
set / of three different 4 x 3 symbols which exist due to the substitution (p qr) = 
(123) cyclically in (17) account for the 144 symmetries exhibited by (14). Itis now 
straightforward to obtain the closed form expression for the 6—/ coefficient : 


eae el Serena ee Scape 
defs y I(6-A BC 2 ea Ss ee 


(F + G— C — D— 2)/2 
...(19) 
(F + G — A — B—2)/2 
Setting one of the numerator parameters of the 4F3 (1) to — 1 (say, D= — 1) and 
replacing the negative parameters A, B and C by — v, — w, and — u, respectively 
and letting * = y and G = x, we get for ther. h. s. of (19): 
{ (x + w)/2 (y + u + w— 1)/2 (x + y+ u — 1)/2 
(x +u +v — 1)/2 (vy + v)I2 Gabi emery 
...(20) 


which is a symmetry of the parametric solution of Brudno and Louck", in the notation 
of Srinivasa Rao et al? for the polynomial zeros of degree |. 


It was shown’ that the 6—/ coefficient can also be represented by aset II of 
four 4F3 (1) sas: 


& b e : 4'--2 ’ a] : 
depp ted 2NT(4‘)(P (1 — BL 1—C,1 — DE’, F,G)]-1 


* 4F3(A;, BC. DEF Gat) acta) 
where, using (15), the numerator and denominator parameters can be shown to be: 
A 


I 


Rag + Rry 7 Reg + 2, B’ Tat Rp, C oe ace Rpg, dD’ = — Rp3 


jh 


I 


Ray — Rp, + 1, F’ = Ry — Ro, + 1, G' = Rey — Rp paras 
...(22) 
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for (pqrs) = (1234) cyclically. The 4 x 3 symbol for this set II of 1 Fs (1) s can be 

written as : 
: Jip pat  o ata De 

fies Eee Bom) EB = Cem |e pf)’ | 

an F’— B’—1 F'—C'—-1 F’'—DpD'—-1 

G— B—1 G-C—-1 G —D'-1| 

eZ A) 

where the positive numerator parameter A’ does not appear in (23) and this 4 x 3 


symbol exhibits only 36 of the 144 symmetries of the 6—/ coefficient which arise due 
to its invariance to 3! column (or all permutations of B’, C’. D’) and 3! row permuta- 


tions (or all permutations of FE’, F’, G’) of [Rie ||. The set II of four IRi lls which 


arise due to the substitution (p q rs) = (1234) cyclically in (22) accounts forthe 144 
symmetries of the 6—j coefficient. As in the case of the set I of 4F3(1)s, in the case 
of this set If of 43 (1) s also we obtain a closed from expression : 


et oa i aac a 


(i Janse a pea) (Fr Gs Creep’ = 2)/2 
(GB = DA)? (RECS? peer) \ 
...(24) 


Setting one of the numerator parameters of the 4F3 (1) to — | (say, D’ = — 1) and 
replacing the negative parameters B’ and C’ by — x and — y, respectively, and letting 
E', F’, and G’ be v, u and w, we get for the r.h.s of (24): 


ee fm ar (xtutv—1)I2 Ytutw— I)/2 
(x+y +u —1)]2 (x + w)/2 (x + y)/2 
...(25) 


which is a symmetry of (20). We have shown elsewhere!2 that the polynomial zeros 
of degree one of the 6—j coefficient are obtained when the parameters in (25) are 
subject to the condition : 


xyz=uvw ...(26) 
or equivalently, ABC = EFG, D= — 1 in the case of (19) and A’ B’ C’ = E' F' 
G’, D’ = — 1 in the case of (24) which is a multiplicative Diophantine equation of 


degree 3 subject to the constraint : 
Fee Sado te UY OW, tl 2y) 


Obviously, polynomial zeros of degree n arise when the sum of the first » + 1 terms of 
the 4F3 (1) occurring in (19) or (24) adds to zero. We could use either (19) or (24) to 
generate the complete set of zeros of degree n. 
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The polynomial zeros of degree 2 of the 3—/ and the 6—/ coefficients in terms 
of these closed form expressions have been studied by Louck!?® ef, a/. using their con- 
nection to Pell’s equation. However their study does not lead to all the polynomial 
zeros of degree 2 of the 6—/ coefficient. Simple algorithms based on the principle of 
factorization of integers have been proposed by Srinivasa Rao and Chiu!4 to obtain 
all the polynomial zeros of degree 2 of the 3—/ and the 6—/ coefficients. 


In conclusion, we have shown in this article the connection between sets of p41 Fp 
(1) s and sets of Regge or Bargmann-Shelepin symbols for the 3—/ andthe 6—/ 
coefficients. This led us to closed form expressions for the polynomial zeros of degree 
n of these coefficients. As 1 increases, the complexity of the constraint equation which 
has to be satisfied by the parameters in the closed from expressions-or, the numerator 
and dedominator parameters of the »41Fp (1) s—increases. At present detailed studies 
have been made only of the polynomial zeros of degree | and 2 of the 3—/ and 6—/ 
coefficients. 
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In this paper the basic equations of linear micropolar elasticity in polar 
coordinates are arranged in the form of vector matrix differential equation 
in the Hankel transform domain. The problem is then converted to an 
algebraic eigenvalue yroblem and solved inthe same domain. Itis seen 
that the results obtained by eigenvalue approach are in full agreement with 
those of other researchers. Further it seems that this approach is more 
elegant and it is believed that this technique has not been applied earlier 
by any researcher to solve the fundamental equations of linear micropolar 
elasticity. 


INTRODUCTION 


In recent years a detailed exposition of the linear theory of micropolar elasticity 
has been given by Kuvchinski and Aero!9, Eringen and Suhabi9, Eringen8, Nowacki!2 
etc.. Nowacki!! has shown that the equations of motion of linear micropolar theory 
in polar coordinates can be decomposed into two mutually independent sets of three 
equations each in case of axisymmetric problems. Dhaliwal and Chowdhury® have 
solved the set (2.6) for the axisymmetric Reissner—Sagoci problem and the solution 
is obtained by the classical method. Dhaliwal’ has solved the set (2.7) for the axisym- 
metric Baussinisq problem. Das et al.2"> have applied recently the eigenvalue approach 
in solving the basic equations of thermoelasticity and extended the approach to 
magneto-thermoelasticity. They have solved the basic equations by representing them 
to single vector matrix differential equation and converting finally to an algebraic eigen- 
value problem. 


In this paper we apply the technique of eigenvalue to solve the axisymmetric 
equations of linear micropolar elasticity. It is believed that none of the previous in- 
vestigators have applied this approach in solving the problems of linear micropolar 
theory. Here the basic axisymmetric eqns. (2.6) and (2.7) are presented in terms of 
single vector-matrix differential equations in sections 3 and 4 respectively. These lead 
to eigenvalue problems (3.8) and (4.4) respectively for the sets (2.6) and (2.7) and these 
are solved for displacements in the Hankel transform domain. The characteristic equa- 
tion of (2.6) gives repeated roots while the set (2.7) gives real distinct roots i. e. 
real distinct eigenvalues. The general solution for distinct roots is obtained by usual 
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procedure (see appendix A), while for repeated roots the solution is obtained follow- 
ing the procedure of Das etal.2 (see Appendix A). Further, solution for the half- 
space is also obtained. It is also observed that the solution obtained by Dhaliwal®”? 
for sets (2.6) and (2.7) are in full agreement with those obtained by eigenvalue approach. 


2. THe Basic EQUATIONS 


The equations of motion and other basic equations for a homogeneous isotropic 
centrosymmetric linear elastic bedy occupying a region R (vide Dhaliwal’), are given 


by 
Oji,j + Pxi = Pitt | (2.1) 


pgt,j + etjk afk + Jyt = J ot 
and the kinematic relations are 


Bis = yt 


Pa pe ae 
\yig = Ujyt + Ek Wk 


the linear constitutive law being 


otj = A yee dey + 2wyeyy + 20 ytAsl 3) 


pig = B3Bxe dty + 2y Baz) + 2e Bras] 


where ctj are the stress tensor ccmporents; ij the ccuple stress tensor components; 
ut the displacement field components; «4 the rotational field components; X% the body 
force components; ¥i the tody ccuple components; y1j the strain tensor components; 
Bi, the curvature twist tensor components; «ijk the unit antisymmetric tensor; [ ] and 
( ) indicate respectively the skew symmetric and symmetric part of a tensor; A, », «, 
B, y,«, are the elastic constants of the micropoljar material; P isthe density; J the 
rotational inertia; and the dot (. ) denotes the derivatives with respect to time. 


Substituting (2.2) and (2.3) in (2.1) a set of six differential equations are obtained 
and these equations are presented in the vector form as: 


. 


(i) (A+ Bw yy.u—(@utadVv xv x u+Z7ay x meee 

(i) (B+ 2y) VV. w~(y +6) VXAXxw + 209 x u— 4aw+J¥= Jo 
...(2.4) 
Here we observe that the material constant « is responsible for a coupling of a displace- 
ment and micro rotation fields. Though these equations are coupled, they are inde- 
pendent in the case, when « = 0. In this case eqn. (2.4); (i) reduces to displacement 
equations of motion of classical elasticity and eqn. (2.4) (ii) describes a hypothetical 


clastic body in which only rotation occurs. When « > eo, the couple-stress theory 
conditions 


w= $VxXu 
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is obtained. 


Now to solve the static problem with no body forces, we take ¥ = Y = u = w 
= 0 and the cylindrical polar coordinates (r, 9, z) is introduced. Equations (2.4) now 
assumes the forms 




















J Ur 2) Otls de 
(4 + a)(v? x hey 2 6 J+ a+ pay 
1 da, 8 wz - 
ee eae 
u 2 Ou 1 de 
(1 + 4)( V2 — “% =f 2 a) +Q+u- a4 2 @ 
Pe 2ir ets Pishaael & 
2 ra ame 0 
> de ] é dw, 
(4 + a) y2uz + (A + p — a) aa + 2a —— wri i > Dia a0. 
= 0 
pad 
(+9 (ve or - 3 a ast )-4eor + @+y-0 





oy a) pisat Ustee 0 Ue }=0 
+2a(— ao Oz 


@ Die Ons sae OX 
( +0(v2o, aor Fee) — tame +B +9 €) dQ 


6 ur ou) = 











(+9 Vor deur +B+y—0 & te A (rug) 
aa | a, (2.5) 
where 
ee tet ae 
vies + ae Pie 7 oo AZ) 


The case in which the vectors of displacement u and rotation w depend only on the 
coordinates r, z and as such the equations (2.5) are decomposed into two mutually in- 
dependent set of equations, viz. 
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de OMe _ 
(a) (1 +0) (vim - 2) +t na) ap 2% Par aaa 





de fe) al \ 
(b) (u +a) y2uz2 + (A +p — a) + 2«—— yp Are) = 0 








w® Cur 6 Uz os 
(c) (r+) (vPaq — 22 ) — dawg + 20 (- ae )= 

















é er J 
(2.6) 
5 Ug OW, Cw, yi : 
(» + x) ( v2 up —  ) 420 Pais ie = 
oy 
(+ 9( v2 oF - or) — anor t+ B+y—9 2 | 
, ou 
— 2aae =(0 (y+ e) V2wo — 4daws + (B+yv—e) or | 
, ' 
+ 2a——(r Ue) = j 
ey af 
where 
l é Cuz 
o> r or (rue + oz 
] rs) rs) Gz 
a | rar ay Oz 
o . & [oe 2 
eer Tees ee Y 


3. SOLUTION OF EQUATIONS (2.6) 


Here only the set (2.6) are considered. The following state of force 


couple stress » are being ascribed to the displacement vector u = 
tion vector w = (0, « a 


-stress o and 
(ur, O, vz) and rota- 


Orr 0 Crz 0) 


c = 0 


Mre 0 
Sep O} p= Mer O * fez 
Szr O azz 3 Hep 0 siGas8) 
where the particular components of stress-tenso 


T have the following forms after using 
the relation (2.3) 


a 
Sry == Qn = + re, 





u 
07? = 2p ro re 


é 
Cis ea 2s are eae 
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eu du Q a 
oe (E+ St) a ( eK te) 2a wy 











6z or 
Oulz our) Oty Guz 








Owe We Ior9 Qe 
= ¥ — - € 
Pre 20) ( Or ) +: ( ap + : ) 
( OM~ Wy ) ( GW re We 
= > = S75 reak =e ae Siu 
Per Y op - hay : 


Cw é r 
Bee = (r — ©) —3*, pee =(¥ +9) 7 (5,2) 








In the system of eqns. (2.6), three mutually independent functions uy, vz and we are 
involved. Multiplying (2.6) by Jo ( r) and (2.6a,c) by Ji (Er) and integrating between 
the limits 0 to oo, we find that the system of partial differential equations (2.6) reduces 
to the following system of ordinary differential equations: 


. ) 
(4 + 9) DP ~ A+ 2p) Bar +p a) ED ie — 24 Dog = 0 | 
Ou a) OD. a, > [(A + 2p) D? = (u + x) Bile + 2x Boe = 0 + 
- | 

2a Diy + 2a € tie + [(y + €) (D2 — &) —4 4] 0? = 0 J 
ra) 


where dy, iz and Og are the Hankel transforms of the functions u;, uz and ©? respec- 
tively and are given by 


« 


(lr, 9) =f (urs o~) & Ja Er) dr 


iz = JS uz &Jo(ér)dr 
0 


and 





oe ee 
D= dz ,Dr= dz” 


In the matrix rotations eqns. (3.3) may be represented as 


pta 0 0 7D2f tw |-[ 0 (A+p—a)& 2a 
0 A+2p O | liz —(Atp— ae 0 0 
0 0 yte Wy — 2a 0 0 


(equation continued on p. 1242) 
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D hae (A+ 2u) & 0 0 te 
liz 0 (u + a) &2 — 208 tie | —Q 
Oe | 0 —2a& (y+e)E2+2a 11 we 
...(3.4) 
we write 
pt « 0 0 plier (A+tp—ay& 2x 
M= O A+ yp 0 — (A+ p—a& 0 0 
0 0 ¥ +« | — 2a 0 0 
7 = [- ly > xX’ aaa D ty 
Uz iz 
af a 
(A + 2p) 5? 0 0 WX Sera By 
a= 0 (anfogh es —2ak itz 
0 — ak (y + €) E2 + 4x We 
Bi = M-1N 
and 
Bo = M71 Pp, A ERY) 
Equation (3.4) with the aid of (3.5) reduces to the form 
Xe se Ba” ab, Bo X. ASO} 
Using block matrices eqn. (3.6) assumes the form 
Cie? % p Bi Be 2 G.7) 
A ii ella ie reas 
where / is a 3 X 3 unit matrix. 
Writing 
By Bo x 
A= [ | ana V= | ] .-.(3.8) 
I Oo X 
eqn. (3.7) assumes the form 
d 
a mee ...(3.9) 


Assume V = Y exp (¢ z) to be a solution of eqn. (3.9). Then we must have 
AY ot LY, (3.10) 
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This gives that ¢ is an eigenvalue of the matrix A and Y the corresponding eigenvector. 
The eigeuvalues of the matrix A are the roots of 


det (A — tJ) = 0. 











<a Gdeu) 
Using simplified notations, eqn. (3.11) may be written explicitly as 
aah aj2 a3 a14 0 ‘ior 
sh 0 0 425 a6 
431 0 “al 0 435 436 0 
ot 0 0 ...(3.12) 
0 —t 0 
0 l 0 th 
where 
2A be =e yl fra A + 2p 
a12 aa &, 413 = eee a14 > Petre e2 
At eH oF ys a lo 7 2a8 
Ly = 28 cat —2a%& — (y + 28 + 4a 
a31 = pri a35 = tent ee Ae » ahs 
eto. hs) 
Simplifying (3.12) and using (3.13) therein we get the characteristic equation as 
16 — (282 + C2) 24 + (E4 4 282 2) 122 — F402 = 0 ...(3.14) 
where 
v2 = E24 m2 
and 


i ap 
~ (w+ a) (y + 6) 


m2 
The roots of (3.14) are 


Jo iN Ae 9 ip os 


The eigenvalues of the matrix A are the root of the equation (3.14). Write ti = &, 
lo = — &, tg = C, and ta = — ¢. The four eigenvectors corresponding to four dis- 
tinct eigenvalues 11, fo, f3, fa of the matrix A are obtained by solving the following 
homogeneous equations 
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f —t @e 43 a4 0 0 } f yi (t) | 
| a21 —f 0 0 a25 426 | ye (t) 
| @31 ae oe! 0 435 436 | | y3 (1) 0 
| 1 0 0 —t 0 0 | | va (t) : 
| I if cs (Oy aKeteam (One 
0 I ad ee ; [ ye (t) j 


(3.15) 


for t = ft, to, tg tq. 


Denote by A; (t), Ag (1), ..., Ag (t) the co-factors of the elements of the first row of 
the coefficient matrix in eqn. (3.15), then 


Ay (t) 1 
Ap (t) | 
y(r)= | 43 (4) : 
Aq (t) | 

| 
As (t) | 

J 


3.16 
Ag (t) 3,16) 


Ce ee ee ee mane, | 


are the solutions of eqn. (3.15) and hence they are eigenvectors corresponding to the 
eigenvalues f1, fg, f3, t4 of the matrix A. By actual calculations we obtain 























= eu ios eerie = 7) 
Aas : ee om = (28 = ES) (A Tete | 
io - os 5°) : 
a= te] SESS py ane |- oe ere 
Ao (0) =) Enea ..(3.17) 


Since tf; and fg are double root of the characteristic equation (3.14) of the matrix A, 
the solution of the differential equation (3.9) is given by (vide Das er al.?), 
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V = Ci Y(ti) exp (t1 z) + C2 d/dz [Y (t) exp (t Z)]tat) 
+ C3 y (te) exp (t22) + C3 d/dz[Y (t) exp (tz)}--t, 


+ C5 Y (ts) exp (tg z) + Ce Y (ta) exp (t4z) ...(3.18) 


where Cj, Co, ..., Cg are arbitrary constants to be determined from boundary con- 
ditions. 


Equation (3.18) can be rewritten as 
V = (Cy + Coz) Y (ty) exp (41 z) + Co Y (t) exp (fz) 
+ (C3 + C4 z) ¥ (ta) exp (to z) + Ca Y (to) exp (foz) 
+ C5 Y (t3) exp (tg z) + Coe Y (ta) exp (t4z) e321 9) 


where dot (. ) represents the differentiation with respect to t. 


For the half space z > 0, equation (3.19) reduces the form 


V = (C3 + Cg z) Y (te Zz) exp (tg z) + C4 Y (to) exp (f2z) 
+ Ce Y (ta) exp (ta 2) ...(3.20) 
where the constants C3, C4 and Cg are to be determined from the boundary conditions. 


Equations (3.20) can be written explicitly as 
a’ (z) = {(C3 + Caz) Ai (to) + Ca Ai (t2)} exp (12) 
+ C6 Ax (14) exp (taz) 


a (z) = {(C3 + Ca z) Ao (ta) + Ca Ao (t2)} exp (t22) 
+ Cg Ag (ta) exp (t4z) 

@, (2) = (C3 + Caz) As (t2) + Ca As (t2)} exp (122) 
+ Cg Az (ta) exp (faz) 


iy (z) = {(C3 + Caz) Aa (to) + Ca Aa (t2)} exp (t2z) 
+ Ce Aa (ta) exp (f4z) 
iz (z) = {(C3 + Cy z) As (fo) + Ca As (t2)} exp (t22) 
+ Ce As (ta) exp (ta z) 
@y (z) = {(C3 + Caz) Ao (t2) + Ca Ag (t2)} exp (f22) 
+ Cg Ag (ta) exp (t4z) ...(3.21) 


1246 R. K. MAHALANABIS AND J. MANNA 


where Cz, C4 and C¢ are arbitraiy constants to be determined from the boundary con- 
ditions. Thus the displacements have been obtained in the transformed domain and 
as such the stresses can also be obtained frcm (3.2) using (3.21) in the transformed 
domain. 


4. SOLUTION OF Equations (2.7) 


Here we are concerned with the set of equations (2.7) in which the displacement 
vector u = (0, ug, 0) and the rotation vector w = (@;, 0, wz). The field of displace- 
ments (0, vu», 0) and rotations («, 0, «z) described by the set of equations (2.7) induces 
the following state of force-stress and couple-stress (vide, Dhaliwal) 




















6s = 0 ore 0 and i err 0 -rz 
Sor 0 Sz 0 Too 0 ..(4 1) 
0 Sze 0 Per O pez 
where 
ou 
ore = (» + a) S* — (H — a) —t - 2s 
a 
Ser = (up — 2) > — (peta) “% + 2 w, 
r r 
a 
Gez = (wp — a) 2 — 2a we 
r) Ug 
Czy = (uw + a) PE + 2% wy 
Wee Peay OW, am 0 Wz 
i Y gp Bee = BX + 2r az? Mee BX + 2y or 
ad OWws ew 0@ 
bre = (y — €) az + (y¥ + €) ae > per = (y + «) a 
‘ 0 wz 
oe Cyoes ig) ce k . ...(4.2) 


Now Hankel transform of the set of equations (2.7) give 


[(y + «) D? — (B + 2y) & — dalor —GB+y—)ED a, 


— 24 Dip = 


(B + y — JED wy + [(B + 2y) D2 — (y + €) E2 — 4 a] oe 
+2a€i% = 0 


—— __—— + 


2a Dw, +28 w, + (u + a) (D2 — £2) i, = 0. 
me) 
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Now as in section 3 equations (4.3) can be written as vector- matrix 


tion form as 
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differential equa- 








d 
mek WY BY (4.4) 
where 
Vea f ow! ; B= ; 0 bie big bi4 0 0 } 
| | | 
wey I | bez O 0 0 bo5 bag | 
| w' | | 
| : bi 0 0 O bss bse | 
teen | | 
| ve | | l 0.6 10m209. 0 0 | 
ae | | 
Were Mees ee (et ye 
lpi 
ae DOR tas Ey re oe 
it uy J Xe J 
...(4.5) 
and 
+y—e, vr, (8 + 2y) &% — 4« 
b > = petnys=. = = ~ = 
12 ae ees Bh ie bia Sta 
= —_ 24 4 —2a 
eee ee oat ee ¢ Ce ae yp a 
21 Be 2y E, bas B+ 2y , dr rare 
= oe —2a€ 
baj = b = ——_, hg = &2 --.(4.6 
31 See hae catia mee E (4.6) 


Assume V = X exp (tf z) be the solution of the equation (4.4). Then we must have 


BX —t X. 


(4.7) 


This gives that ¢ is an eigenvalue of the matrix Band X the corresponding eigenvec- 


tors. 
det (B—1t/) = 0 

That is, 
— f bie big bi4 
bo1 ay 0 0 
bs 0 —t 0 
| 0 0 —t 

] 

0 0 0 


where by2, bi3, etc. are given by (4.6). 


The eigenvalues for the matrix B are the roots of 


.-.(4.8) 
0 0 
bes 26 
b35 36 
0 0 = 0 (4.9) 
—-t 0 
0 -t 
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Now simplifying (4.9) and using (4.6) therein we obtain the characteristic equation 
as 


16 — 74 (# + art +2 Jee (ea; + B25 + ai ab ) 
pL} ’ 
— BAAS =0 ---(4.10) 
where 


A= +Ki, B= e+ K? 


I 


4a 2 4ap 
Ko = - | Oe Ed we 
1 B =. 2y- ° 2 (u + a) (y + €) 


whose roots are €,—£, Ai, — Ai, Ao, — Az, which are the distinct eigenvalues of the 
matrix B. The corresponding eigenvectors are obtained by solving the following homo- 
geneous equation. 


( —t bis big big 0 0 : [ SK ED, 1 
| bey —t 0 O bes bog | | x2 (t) | 
| 63, O —t O h35 bag | | xs(t) | 
[pees | =O. ...(4.11) 
| ] 0 QO -2 0 0 | fe fear (1) | 
| 0 l —T 0 | , x5 (t) | 
RD kl P0900. pe Beery | 
fort = 4, i= 1, 2,..., 6 wherety = &, to = — & ts = Ay, t4 = — Aj, ts = Dz, 


16 = — Ao. 
Denote by A: (t), i = 1, 2, ..3; 6, the co-factors of the elements of the first row of 
the coefficient matrix in eqn. (4.11), 
Then 

f By (t) 
By (t) 
Bs (t) 
Bq (t) 
Bs (t) 
\ Be (t) 


are the solutions of the eqn. (4.11) and hen 
the eigenvalues 4, i = iva 


| 
X (t) = | 
| 
| 
| 
i 


— 


ce they are eigenvectors corresponding to 
.» 6, of the matrix B. By actual calculation we obtain 


A= 84a) 


(equation continued on p. 1249) 
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—1® (eno tee 
(8 + 2y) (H+ 9) j 





fia (B + y — «) aaa - 
Bo(t) = — i: a (u + a) (B+y — a€2) 
(8 + 2y) ‘aaa ret Taal CS a) } 
B3 (t) = — 14 ee a ee Oe 
ee (un + «)(B + 2y) 
B(y=n- 2] Etre | 
Y 


tye ahoneh) A Rawk caer 8 
(4 + 4) (8 + 27) 





Bs (t) = 13 Bty—+_ {fra dot 
B + 2y Speers Bie GES tice \ 
Cal ea sae ~ p(B tv =) — rely FOE + bet 

is (B+ 2) +4) t 


Since ft), fo, ..., fg are all distinct roots of the characteristic equation (4.10) of the 
matrix B, the general solution of the differential equation (4.4) is given by (vide 
appendix A) 


V = £, X(t) exp (t1 z) + Eo X (ta) exp (fe z) + E3 X (ts) exp (ts z) 
+ Ey, X (ta) exp (ta z) + Es X (ts) exp (t5 2) 
+ Eg X (te) exp (t6 Z) 
where 
To ae fae eg rly a a 
are arbitrary constants 


1. e€. 


6 
V= > Ei X(t) exp (tt 2) 
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PINAY 


APPENDIX A 


Consider the differential equation 


dy 
a / = Ay .. (Al) 
where ¥ is ann vector and A is ann X n real constant matrix. 


If Ay Ag, ..., An are distinct eigenvalues of the matrix A and x1, x2 ... xn be the 
corresponding eigenvectors of A, then the general solution of (A1) is given by 


v(t) = Cy Xy ey! + 2... + Cn Xn eAnt 
where Cj, ... Cn, are arbitrary constants. 


If Ay is an eigenvalue of A of multiplicity 2 and all other eigenvalues A3; ... An 
are of multiplicity one and YX, X3,... ¥n the corresponding eigenvector of A, then the 
general solution of (A!) is given by 


f; es 
v(t) = Cy Xy-e444 + Co ae (Xz e44") + Cz X3 erst + ... 


+ Cn Xn ernt 
where €y, Co,..., Cn are all arbitrary constants. 


For details, vide Das et ai.2. 
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The Mellin transform of the gravity effect of a buried 2— D horizontal circu- 
lar cylinder with density contrast varying linearly with depth is derived and 
analysed to extract the body parameters. The characteristic features of 
the Mellin transform of the gravity effect resemble gamma function. The 
validity of the method is tested on simulated models. The stability of the 
Mellin transform is studied by incorporating random noise in the gravity 
effect at different levels and subsequently the error estimation of the inter- 
preted values is discussed. 


INTRODUCTION 


In general, the interpretation of the geophysical anomalies are carried out by 
a ssuming certain geometrical shapes with uniform physical properties. Very frequently, 
we come across some practical cases wherein for example the density contrast increases 
with the increase of depth which is evident from seismic studies. 


Therefore, it would be meaningful to interpret such geophysical anomalies with 
non uniform density4. Herein the analysis of gravity anomalies due to a horizontal 
circular cylinder with the density contrast varying linearly with depth is presented us- 
ing the Mellin transform. Such transformation of gravity (or magnetic) anomalies 
paves the way for simplified analysis of the complex potential field!. The procedure 
is illustrated with three sets of theoretical models. The stability of the Mellin transform 
is studied by incorporating random noise in the gravity effect at different levels and 
subsequently the error estimation of the interpreted values is discussed. 


MELLIN TRANSFORM OF THE GRAVITY EFFECT OF THE HorRIZONTAL CIRCULAR CYLINDER 


A buried horizontal circular cylinder extending infinitely along the Y-direction, 
with its normal section parallel to the ¥ — Z plane is considered. The origin of the 
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coordinate system is taken on the ground surface such that the Z-axis coincides with 
the diameter (Fig. la). Let the density contrast at the apex of the cylinder as P and 
the rate of change of density contrast varying linearly with depth be a. In this case, 
the gravity effect of the cylinder is given by Radhakrishan Murthy. 


Z Zi: Se Bye 
EQ) = 4 Yat gs — 8 ya gape 2{)) 
Where 
A = 2n GR2 (p + aR) ...(1a) 
GaR4 
SOS mee areas (1b) 


Z is the depth to the centre of the cylinder, R the radius of the cylinder andG the 
universal gravitational constant. 


The Mellin transform of the gravity effect given by eqn. (1) is written as Sneddon® 


M(s)= f xe 
(s) J x81 g (x) dx Az) 


300 


HORIZONTAL CYLINDER MODELS 
(With variable density ) 


g(x) 


200 


IN ARBITRARY UNITS 


100 





Se eel) > 


Fic. 1(a) Cross section of the buried horizontal circular cylinder. 


(b) Computed gravitational effect of Norizontal circular cylinder 
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where s is a real positive integer or fractional number. 


Substituting for g(x) and integrating eqn. (2) with respect to x, we get the Mellin 
transform of g (x) as: 


M (s) = [A Z#-1 T(s/2) T ((2 — s)}2) — BZs-2 
x {P (s/2) T(4 — s)/2) — P(s + 2)/2) 1(2 — s)/2)}] -..B) 
(Oia ge <7 2) 


ANALYSIS 
For a set of arbitrary values of s, i,e., for s = 1/4; s = 1/2; s = 3/4; 5=1 
and s = 5/4 equation (3) is written as : 


M (1/4) = A Z73/4 1(1/8) 17/8) — B Z-7/4 [ PU/8) T(15/8) 


— 19/8) T7/8)] ...(4a) 
M (1/2) = A Z-1/2 T(1/4) 1(3/4) — BZ 3/2 [ Pd/4) 17/4) 
— T(5/4) P(3/4)] ...(4b) 
M (3/4) = A Z 1/4 1(3/8) 1(5/8) — B Z5/4 [ 1(3/8) (13/8) 
— (11/8) 1(5/8)] ... (4c) 
A= M(\)/a ...(4d) 
M (5/4) = A Z1/4 1(5/8) 13/8) — BZ-3!4 [ 1(5/8) M(11/8) 
— 1(13/8) 1(3/8)). ...(4e) 
From eqns. (40) and (4e), the value of Z is evaluated as: 
M (3/4) Z1/4 = U — (BIZ) V ...(5) 
and 
M (5/4) Z1/4 = U + (B/Z) V tO) 
lec, 
P(z142 —U Z1/4+Q=0 matt) 
where 
P = M (3/4)/2, 
Q = M (5/4)|2 
U = M (1) (3/4) 1(5/8)/7; 
V = (3/8) T(i2/8) — T1178) T'(5/8). 
Hence 


U + (U2 — 4PQ)i/2 \4. ...(8) 
ae, 
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Since A and Z are known, B can be evaluated as : 


-1/4 — 1/4 
Beez | Oe ey 1, (9) 
By eliminating ‘a’ from equations (la) and (Ib) a cubic equation in ‘R’ is obtained as : 
R3 — (A/2 P 7G) R + (4B/2 P nG) = 0. .. (10) 

Apply ing the well known Cardon’s metked, Ris evaluated, and subsequently a 


is calculated as: 


a= 2B/Ré. ms ait 


DiscRETE MELLIN TRANSFORM 
Since the gravity data is collected at discrete intervals in the real field situation, 
the numerical computation of the Mellin transform is carried out by formulating the 
discrete Mellin transform as? 






HORIZONTAL CIRCULAR CYLINDER 
(With variable density ) 





Theoretical Metin Transform 





"4 7) Wa 1 Yq ¥ % 2 


—— S 
Fig. 2. Continuous Mellin transform of the 


gravity effect of i : : ; 
density variation. y of horizontal circular cylinder with 
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N-1 
deemed Eden g(n.& x) (n, Ax)!.As—1l Ax vata) 


(Oe is As < 2) 


where As = 1/4, N is the total number of points and Ax is the sample interval. 


SYNTHETIC EXAMPLES 


The procedure detailed in the text is illustrated with three theoretical models 
(Table I). The theoretical Mellin transform of the gravity effect of three models are 
computed using eqn (3) and shown in Fig 2. 


The gravity effect due to horizontal circular cylinder with variable density for 
three models are computed using eqn. (1). Since the gravity effect due to horizontal 
circular cylinder is symmetric, only the positive side of the anomaly is shown in Fig. 1b. 
The discrete Mellin transform of the gravity effect of the cylindrical models are com- 
puted using eqn. (12) and shown in Fig. 3. It may be observed that the discrete Mellin 


14 HORIZONTAL CIRCULAR CYLINDER 
(With variable density) 


10 Discrete Mellin Transform 


M(n-AS) ——> 









7 7 
0 \, VY, ¥, 1 S/, % Ms 2 


(VAS) — 


Fig. 3. Discrete Mellin transform of the gravity effect of horizontal circular cylinder with den- 


sity variation. 
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TABLE | Theoretical Examples 
R* h* a* 
f Assumed values 0.50 ON ba 1.00 
Model I 
Evaluated values 0,56 0.70 0,92 
Assumed values 1.00 2.00 1.50 
Model II 
Evaluated values 0.95 1.88 1.39 
Assumed values 1.50 3.50 TS 
Model III 
Evaluated valves 1.45 2.95 1.80 
(* in arbitrary units) 
| 
Model —I 
10 cates 
0-9 _ 
0:8 —e 10% 
x X 20% 
0-7 o---0 30% 
b= 8) 40% 
0-6 
0-5 
a. 1-1 
c Model — I 
@ 
= ie mete 
= 09p 
= SRE 
CO 08 
= 
o 07 
= 
© 0:6 
o 05 
© 
1-2 . 
ie} Model — It 
1-0 
0-9 
0-8 
0-7 
0-6 ; ; 
100 200 300 400 500 600 700 
Number of samples ——»> 


Fic. 4. Correlation 


coefficient 


versus number of discrete gravity samples, 


_ 
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Model —I 


——2 -10'% 
Joserevee) X 20% 
o--9 30% 
&--b 40% 


Model —II 


<————_ Standard deviation 


Model—II 


—_— 
= 
= 





108 200 300 400 500 600 700 800 
Number of samples ——» 


Fic. 5. Standard deviation versus number of discrete gravity samples. 


transform of the gravity effect of the models (Fig. 3) are similarto the theoretical 
Mellin transform (Fig. 2) and they resemble gamma function curves. The parameters 
are evaluated from the computed discrete Mellin transform of the gravity models 
using eqns. (4d), (8) (10) and (11) and Fig. (3) and presented in Table 1. It may be 
noticed that the evaluated parameters reasonably agree with the assumed values. 
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12 Model-—I os 10% 
X----% 20 % 
me) o--0 30% 





Model —II 


variation —— 


error 


~<--—— Normalised 
°o 
a 


Model —-II 





300 400 «= 500S—s«&6OD.-—«<0-S~SOD 
Number of samples —~e 


Fic. 6. Normalized error variation versus number of discrete gravity samples. 


STABILITY ANALYSIS 


Here for all the three gravity models the number of discrete samples of the gra- 
vity effect are considered at six different levels ranging from 100 to 600 at an interval 
of 100. In each case the random noise with different percentages i.e., 10% 20%, 30% 
and 40% are added to the discrete gravity effect of the horizontal circular cylinder 
with variable density. The discrete Mellin transfom of the noise contaminated gravity 
effect is computed. Using computed theoretical Mellin transform of the gravity effect 
and the discrete Mellin transform of the noise contaminated gravity effect forall the 
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models for different levels of number of samples and with different levels of random 


noise, the correlation coefficient, standard deviation and the error variation (normalized) 
are computed and shown in Figs. 4, 5 and 6. 


It is observed from the Fig. (4) (i.e.) correlation coefficient versus number of 
samples with different levels of noise for the gravity models, the correlation coefficient 
is tending to | for discrete gravity samples > 400. Also the standard deviation and 
error variation (normalized) are saturating for discrete gravity samples > 400. 


The error percentage in evaluated parameters of the cylinder is relatively high 
(> 30%) for the discrete gravity samples < 400 for different levels of random noise. 
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